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Summary--The absorption of SO, by water and by aqueous solutions of HCl, NaHSO, and 
NaCl using relatively long wetted-wall columns (12 to 22 cm) with no gas phase resistance has 
been investigated. The results of the investigation show (1) that in the presence of small con- 
centrations of a surface active material (here Teepol, an alkyl sulphate wetting agent) the rate 
of absorption is very well predicted by the penetration theory, (2) that the primary function 
of the wetting agent is to prevent the formation of ripples in a falling film of the heights used 
and (3) that the hydrolysis reaction of the system SO,-H,0 is very rapid relative to the contact 
times involved, so that absorption of SO, by water may be considered as physical absorption. 

The data obtainable with this type of absorber may be used to calculate a rather reliable 
value of the product C* 4/D in the case of physical absorption. If this product is known, the 
absorber may be used to study reaction rates in chemical absorption. 


Résumé—On a étudié l'absorption de SO, pure dans l'eau et dans des solutions aqueuses de 
HCl, NaHSO, et NaCl avee une colonne & parois mouillées, de 12 & 22 cm. de longueur. Les 
expériences démontrent (1) qu’en présence de faibles concentrations d'un détergent (Teepol) la 
vitesse d’absorption peut étre prédite avec la theorie simple de diffusion non-stationnaire, (2) que 
la fonction primaire du détergent est de supprimer les ondulations a la surface du liquide, et 
(3) que la réaction d"hydrolyse du systéme SO,-H,O est trés rapide auprés les temps de contact 
utilisés (entre 0,1 et 1 seconde). 

Pour l'absorption physique, des expériences dans un absorbeur comme utilisé ici peuvent 
fournir une valeur assez sare du produit C* ,/ D. Au cas oi cette valeur soit connue, l'absorbeur 
peut étre utilisé pour l'étude de vitesses de réaction chez l'absorption avec réaction chimique. 


1. INTRODUCTION 


end- and begin-effects, inexact methods of analysis 


The penetration theory of absorption in the 
absence of gas-phase resistance, first proposed by 
Hicsie [1] in 1935, has long been recognized as 
being intuitively more reasonable than the 
stagnant liquid film concept which is used so 
widely in practical design work. While many 
investigators [2, 8, 4] have used the penetration 
theory to correlate their data, however, attempts 
to verify the theory quantitatively have not been 
highly successful. Reasons for this have been the 
appearance of ripples on the liquid films, unknown 


and various other complications. As a result, the 
precision of the data has generally been low. The 
present investigation was undertaken in an 
attempt to verify the penetration theory quanti- 
tatively for a very simple flow situation. 

As applied to the case of a falling liquid film, 
the basic assumptions of the penetration theory 
are as follows : Gas is absorbed by an element of 
area on the surface of the film as if the surface 
element were saturated during the first instant 
of exposure, and as if the liquid behind the element 
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were an infinitely deep stagnant layer during the 
time, ¢,, in which the element moves from the top 
of the column to the bottom. Thus, the velocity 
of that part of the film through which the solute 
diffuses during t, is considered to be constant, 
and the diffusion in the direction of flow is 
assumed to be negligible. These assumptions are 
expressed in the following familiar differential 
equation and its boundary conditions 
C=C, forz >0,t=0 
D—=—, C=C* ,2=0,t>0 (1) 
C=C, ,2=0,t< 0 


where D = diffusivity of the solute in the liquid, 


em? /sec, 
C = concentration of diffusing solute, 
g/cm’, 
C* = concentration of solute at saturation, 
g/cm’, 
C, = initial concentration of solute in bulk 
of liquid, g/cm', 
a2 = coordinate in the direction of diffusion 
em, 
t = time, sec. 


The rate of absorption, ®, (g/sec), of an 


element of area at time ¢t is then 


d®,, 


and the average rate of absorptoin per unit area, 
©”... (g/em*sec), for the whole column can be 
shown to be 


» D 


=2(C* — (2) 


D. 
mt, 
where ¢, is the total contact time of an element 
of area. 

For the case of a film formed on the outside 
of a cylindrical tube, the area A, of the film 
surface is 


A = 2n(r + dh, 


where h, = the effective height of the column, 
r = the radius of the tube, 
8 = the thickness of the film. 


(8) 


The thickness of the film is determined by the 
density and viscosity of the liquid and by the 
flow rate. For laminar flow it is given by 


§ = ¢ 4 


where » = kinematic viscosity, cm*/sec, 
I’, = volumetric rate of flow per em, em® - 
cm sec, 
g = acceleration due to gravity. 981 cm,- 
sec?, 


The velocity of the surface of the film is 
similarly determined by the nature of the liquid 
and the flow conditions and for laminar flow is 


The flow has been found to be essentially 
laminar [5] for values of Re below 1200 when Re 
is defined as 

v 
The time during which an element of area on the 
surface of the film is in contact with the gas as 
it moves from the top of the column to the bottom 
is 
r, 


(6) 


Thus, the average rate of absorption per unit area 
for the column becomes 


3 DI, 


Dancxwerts [6] has provided a means of 
estimating the conditions under which the 
assumptions of an infinitely thick layer moving 
at a constant velocity are no longer approximately 
valid. From a graphical solution of the differential 
equation for diffusion in a finite film moving with 
a parabolic velocity profile he found that the 
effect of the wall and the velocity gradient did 
not change the rate of absorption by more than 
5%, so long as the value of the dimensionless 


INLET GAS SATURATOR 
ROTAMETER 
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group = is less than 0-1, 
value of I. used in this set of experiments the 
corresponding value of h should thus be less than 
200 em, a condition which is met with quite a 


wide margin of safety in these experiments. 


‘or the smallest 


2. Description oF EQUIPMENT 
A schematic diagram of the equipment is shown 


in Fig. 1. Water was throttled through a valve 
from a pressure of two atmospheres through a 
rotameter and into the absorber. The purpose 
of having the water supplied under pressure was 
to prevent small changes of pressure in the 
absorber from affecting the liquid flow rate. 


OUTLET CAS 
ROTAMETER 


T 
M 
SO, cas "| 


7 = THERMOMETER 
M = MANOMETER 


Fic. 1. Sketch of absorption apparatus. 


INLET WATER 
ROTAMETER 


WATER 


The water flowed into the absorber through a 
vertical tube made of stainless steel with an 
outer diameter of 15mm. The surface of the 
tube was roughened slightly by sand blasting to 
facilitate wetting by water, and the tube wall was 
slightly rounded at the top. The water entered at 
the bottom of the tube and flowed out over the 
top and down the sides as a falling film. The 
liquid level in the receiver at the bottom of the 
film was kept constant, the height of the column 


"© The substitution of eqs. (4) and (5) in this group 
ah, # 
in other words, 
residence time 


jiffusion time must be less than 0-4, 


being varied by moving the tube rather than by 
changing the liquid level. 

SO, from a cylinder flowed first through a large 
rotameter and was then saturated with water 
before being admitted into the absorber. The 
incoming flow rate of the SO, was always greater 
than the rate at which the gas was being absorbed. 
and the excess passed out through a second, much 
smaller rotameter in the outlet line. The flow 
through the smaller rotameter was about 5% of 
the incoming flow. By thus maintaining a small 
steady stream of gas through the absorber, the 
accumulation of insoluble gaseous impurites at 
the liquid interface was prevented. The tem- 
peratures of the incoming gas stream, of the water 
in the saturator, and of the incoming and out- 
going liquid streams, and the pressure in the 
absorber and before the gas rotameter were 
measured for each experimental point. Small 
adjustments of the data were then applied to 
obtain the rate of absorption at 20°C and a 
partial pressure of SO, of 750mm Hg. These 
adjustments will be discussed in more detail in 
another section. 

Ordinary distilled water was used. Preliminary 
experiments indicated that the amount of air 
dissolved in the water as the result of the pressur- 
ization of the feed vessel had no measurable effect 
on the gas absorption. The SO, used was 99-9%, 
pure, the impurities being H,O and SO,. The 
solutions of sodium bisulphite were made from 
the purest commercially available supply of that 
salt. The normality of the solutions was deter- 
mined by acidimetric titration. The sodium 
chloride solutions were made up by weight. 


38. TREATMENT OF Data 


For ease of comparison, the data obtained were 
all adjusted to a constant value of the surface 
temperature of the water film and a constant 
value of the partial pressure of SO, within the 
absorber. The average temperature of the surface 
of the film was assumed to be 1°C above the 
average of the inlet and outlet temperatures (see 
Section 5). The rate of absorption at the average 
temperature of the experiment (16°, 20°, or 24°C) 
was then calculated by assuming the thermal 
coefficient of absorption to be — 1-6%/°C, a value 
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which is based on the penetration theory (see 
Section 4, A) and which was verified by experi- 
ment. Thus, the rate of absorption observed, for 
example, at an average surface temperature of 
20-5° was increased by 0-8°%, to obtain the value 
for 20°. The adjustments made in this manner 
rarely exceeded 2°, and were usually less than 
1%. 

The partial pressure of the SO, in the absorption 
vessel was computed by subtracting the vapour 
pressure of water at the temperature of the water 
in the saturator from the absolute pressure in 
the vessel. 

The concentration of SO, at the surface of the 
film, and hence the rate of absorption of SO, in 
the liquid, was assumed to be directly propor- 
tional to the partial pressure of SO, in the vessel. 
Thus, the rate of absorption observed, for 
example, at a partial pressure of SO, of 755 mm 
Hg was decreased by 0-67%, to obtain the rate 
at 750 mmHg. These adjustments too, were 
usually less than 1%. 


4 RESULTS AND COoNCLUSIONS*® 
A. Agreement of the data with the penetration theory 


If, for a certain gas-liquid system, one has 
dependable values for all of the parameters of 
equation (7), one may test the penetration theory 
quantitatively for that system. Lacking one of 
the parameters, one may use the absorption data 
of the system as a means of determining its value 
by assuming the theory to hold. In the present 
case, a combination of these two procedures was 
employed. 

The terms in equation (7) subject to the most 
uncertainty for the SO,-water system are, of 
course, C* and D. Values of the former found in 
the literature disagree by about 7°, those found 
for the latter disagree by nearly 50%. Since it 
was thus necessary to make an arbitrary choice 
of the parameters used, we chose the following 
procedure. Absorption data were taken with the 
flow rate, column height, and temperature as 
independent variables. A value of C* was selected 
from the literature, and a value of D was then 


* No numerical results have been tabulated in this paper 
but they can be requested from the authors. 


calculated from the data for one height and 
temperature. Since the temperature coefficients 
of D and C* are fairly well known, it was then 
possible to check the penetration theory by the 
variation of the absorption with temperature and 
column height. 

The data on the solubility and dissociation of 
SO, in water were taken entirely from the work 
of Jounstone and Leppia [7]. Their results 
are presented in terms of an ionization constant 
and a Henry's law constant which are given as 
functions of the temperature. Their values for 
these constants and the total solubility differ from 
those given by other authors. They are used 
here, however, because they are consistent with 
each other and can easily be used in calculating 
the effect of temperature on the rate of absorption. 

In Fig. 2 are shown the data for the rate of 
absorption of SO, by a solution of 0-0050 wt.%, 
surface active material (Teepol) in water for 
columns of different heights. The effective height, 
h,, of the column was obtained by measuring the 
actual height, h, and adjusting for the entry effect 
in the following manner: The area of the liquid 
surface on the top of the tube is more than xr* 


=220m 


=l4em 
= 
—— EQUATION (7) 


0.0500 WT %e SAM 
te 

| 


T 


>>? > 


° 


Fic. 2. Rates of absorption for different heights of column. 


and less than 27r°, i.e., the surface area is between 
that of a circle and a hemisphere. It was taken 
arbitrarily to be 3/2zr?, which can be expressed 
in terms of an increase in the height, Ah, of 
jr. The liquid in the receptor, being relatively 
stagnant, was assumed to absorb SO, at a 


negligible rate. 
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The solid lines in Fig. 2 were calculated from 
equation (7). For C* at 20°C a value of 0-0995 g 
SO,/cem® (data of Jounstone and Leprta [7)) 
was taken. It was found that a value of D at 
20°C of 1-46 x 10°° cm?/sec. fits all the experi- 
mental data. This value is in the same range as 
those found by other investigators. Ko.rnorr 
and [8] reported 1-78 x 10°, PeEaceMAN 
[2] found 1-55 x 10°, Wuirney and Vivian [9], 
as will be mentioned later, based their calculations 
on a value of 1-66 x 10° and Groornuts and 
Kramers [10] found 1-4 x 10% cm? sec.* The 
above values have been interpolated, where 
necessary, to 20°C by assuming Du,/T to be a 
constant. 

The scatter of the data in Fig. 2 is remarkably 
low for absorption work. The straight portions 
of curves drawn through the experimental points 
in Fig. 2 would lie within 4% of the respective 
theoretical lines. Since a deviation of the experi- 
mental points from the theoretical curves of only 
2%, would be significant, it is felt that the agree- 
ment of the two over a range of column heights 
of slightly less than a factor two provides a 
sufficiently rigorous test of the theory for the 
purposes of this paper. 

Furthermore, as can be seen in Fig. 3, equation 
(7) equally well predicts the change in absorption 
with temperatures, caused by variations in C*, D 
and pu. 

The departure of the curves of Figs. 2 and 3 
from linearity at the higher flow rates is believed 
to be due primarily to the increasingly greater 


aoc 4 
‘ 
— EQUATION (7) 
+—+1_0-0500 WT. %e SAM 
TT 
+ 
| 
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Fic. 8. Rates of absorption at different temperatures, 


effect of the entrance conditions. This effect is 
relatively more important for the shorter columns, 
as would be expected. The deviation of the 
points from the theoretical line becomes import- 
ant at increasingly lower flow rates as the height 
of the column is decreased. The minimum height 
of the type of column used in these experiments 
was thus determined by these entry effects. 


B. The effect of the presence of Teepol 

Many observers studying absorption and desorp- 
tion with wetted-wall columns have reported that 
at certain flow rates and column heights the rates 
of absorption measured are much higher than 
those predicted by the penetration theory. This 
effect has generally been attributed to ripples on 
the film surface. These ripples are always 
observable when the rate of absorption is ab- 
normally high, and are not found when the rate 
of absorption is normal. and Picroxp 
[3], as well as others, found that the addition of 
a small amount of an alkyl sulphonate (“* Petro- 
wet ’’) eliminated the ripples and decreased the 
absorption. They also found, however, that as 
the concentration of this wetting agent was 
increased the rate of absorption went through a 
minimum, showing an optimum concentration for 
the elimination of ripples. ZeRNovsKAYA and 
Be.Lopoisku [11] found that the presence of a 
minimum depended on the type of surface active 
material used. 

The effect of the ripples on the gas absorption 
depends a great deal on the height of the column 
being used. At low flow rates the ripples appear 
quite close to the top of the column, and the 
amplitude of the ripples is relatively high. As 
the flow rate is increased, the ripples begin to 
form at points increasingly farther from the 
entrance, and the amplitude of the ripples 
decreases. A column of any height could pre- 
sumably be ripple-free at sufficiently high flow 
rates if it were not for the onset of turbulence in 
the film. 

The addition of a surface active material to 
the water lowers the point on the column at 


* The authors’ value also agrees very well with the 


semi-empirical correlation of Scnerper (Ind. Eng. Chem, 
1954 46 2007) and WILKE (Chem. Eng. Prog. 1949 45 218), 
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which the ripples begin to form. The effect of 
the wetting agent used, an alkyl sulphate marketed 
under the trade name of Teepol, on the absorption 
is shown in Fig. 4. The hump in the absorption 
curve of pure water is gradually ironed out as the 
Teepol concentration is increased. At the highest 
flow rates, where even the pure water does not 
exhibit ripples, it is seen that the absorption by 
the water is nearly the same as that of the 


Teepol solutions. 


a 
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Fic. 4. Effect of surface active material on the rate of 
absorption of SO, in water. 


C. The hydrolysis of SO, in aqueous solution 

According to Rornu [12] in a saturated solution 
of SO, in water the SO, is present in three forms : 
unhydrated, hydrated and ionized. JounsTone 
and Lepp.a [7] reported that the concentration 
of unionized SO, in aqueous solution varied with 
the partial pressure of SO, in accordance with 
Henry’s law, and gave values for the Henry's 
law constant which were based on their measure- 
ments of the solubility and the values of the 
ionization constant which they had chosen. As 
the result of measurements made with a packed 
column, Wuitney and Vivian [9] suggested that 
the rate of the ionization reaction might be slow 
in relation to the contact times involved in their 
column. This assumption can be checked by 
comparing the rate of absorption in water with 
the rate of absorption in a solution which sup- 
pressed the ionization reaction. One should find 


a distinctly lower rate of absorption by an acid, 
for instance, if the ionization reaction were fast, 
and very little difference in the rates of absorption 
if the ionization reaction were slow. 

A saturated solution of SO, at normal pressure 
and temperature is about 1-6 M in SO, and about 
0-18 N in hydrogen ion, i.e., about 11% ionized. 
Because the saturated solution is so concentrated, 
it is necessary to use acid or bisulphite solutions 
of high ionic strength in order to suppress the 
ionization by a measurable amount. Therefore, 
absorption measurements were carried out with 
3 N HCl, and 1 N and 0-2 N NaHSO,,. In addition 
solutions of 1 N and 0-2 N NaCl were used in 
order to provide, for comparison, an idea of the 
effect of ionic strength alone. It is assumed that 
NaCl forms no complexes with SO,. All the 
solutions contained 0-0050 wt.°, Teepol. 

The results of these tests are shown in Fig. 5. 
The lines on the figure represent equation (7) for 
water and the sulphite solutions. To make the 
calculations, the concentration of unionized SO, 
in the sulphite solutions was assumed to be the 
same, when expressed in grams per thousand 
grams of water, as it is in pure water, The partial 
specific volume of SO, was assumed to be the same 
in the sulphite solutions as in water, and the term 
Dy was assumed to be a constant. A summation 


T 
TTT + 4 
TTT - 
© O-OSOO WT % SAM. NO ACID 
o - NaCl 


4=18Ocm 


Fic. 5. Comparison of the rates of absorption of SO, in 
water and in solutions of NaHSO,, NaCl and HCl. 
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% of the pertinent physical and chemical properties 

x of the solutions, measured and calculated, is given 
in Table 1. 

— 5 a eer The rate of absorption by the 1 N bisulphite 

— solution is retarded even more than is predicted 

or & on the basis of reduced ionization and a smaller 

le diffusion coefficient. This is not surprising, since 

one would expect a certain salting-out effect by 

the NaHSO, on the unionized SO,, and the pre- 

~~ = dicted change in the diffusivity coefficient is little 

= e more than a guess. The magnitude of these two 

= effects may be estimated by comparison to the 

absorption by NaC! solutions. As is seen in Fig. 5 

2 0-2 N NaCl absorbs at nearly as great a rate as 

ae Se z : water, while 1 N NaCl absorbs at a rate which is 

»§ GZ FF The rate of absorption by 8 N HCI was found 

RE 3 g\~ = %, to be the same as that of water. This is surprising, 

55 since one would expect both a lower solubility 

of SO, (due to retarded ionization) and a lower 

~ x diffusivity (because of a higher viscosity) in the 

ses 25 acid. Rough solubility measurements indicated, 

if cae en however, that the equilibrium concentration of 


SO, in 3N HCI is about 15% greater than that 
p= of water.. It thus appears that the agreement in 
the rates of absorption is simply a coincidence. 
The evidence given above, in favour of a fast 
hydrolysis reaction, is persuasive but not con- 
clusive. The data could all be interpreted as 
being in accord with the concept of a slow 
hydrolysis reaction if it were assumed that the 
diffusion coefficients in the bisulphite solutions 


0-088 


C*,s0, 
(cale.) 
7 
60-0885 
0-0865 


0-0108 
737566 
471507 


Pertinent physical and chemical properties of solutions used. 


were reduced much more than the change in 

< i l while drastic, can not be dismissed, a priori, as 

Additional data, however, obtained by a 

eco e¢ completely separate method, appear to settle the 


question beyond doubt. An apparatus was con- 
structed which permitted the rapid mixing of a 
stream of water and a stream of liquid SO,, and 


2 the subsequent measurement of the electrical 

zz conductivity of the solution obtained. In the 

i css $ 3 steady-state flow obtaining, complete mixing was 
5 ZZ ~ Zz 7. produced in a time interval of 0-1 second, and 
the conductivity was measured for times between 
Sonacer 0-1 and 1-0 second. Using the conductivity data 


of Morcan and Maass [11], it was found that 


| 
53 


S. Lynn, J. R. StraatTemerer and H. Kramers 


even at temperatures as low as 4°C, the hydrolysis 
reaction was at least 90°%, complete within the 
shortest time interval which could be measured. 

Wuirney and Vivian made their assumption 
of a slow hydrolysis reaction in order to correlate 
their absorption data with that of other investi- 
gators. Their data can be correlated just as easily 
by assuming a lower value of the diffusion co- 
efficient than they used. As mentioned above, 
the variations in the values of D found in the 
literature are sufficiently wide to give one a great 
deal of freedom in making a choice of a value to use. 

In the light of the data and arguments presented 
above, it appears that the rate of the hydrolysis 
reaction of SO, in water is fast relative to the 
diffusion process. Thus, one may assume that 
the surface of a water film is instantaneously 
saturated at the equilibrium concentration upon 
exposure to SO,. It follows that the absorption 
of SO, in water may be treated as physical 
absorption. 


5. Discussion oF Errors 


The random error, as can be seen from the plots 
of the experimental data, is quite small, about 
+1%. It arises from uncertainties in reading 
the rotameters and thermometers and in deter- 
mining the height of the column. 

The systematic errors include both the un- 
certainties in the calibration of the equipment 
and also the uncertainties in the assumptions 
involved in the estimation of the surface on the 
top of the column, of the changes in solubility 
due to temperature and pressure (i.e. of the 
corrections made to present the data for a standard 
temperature and pressure), and of the surface 
temperature of the film. 

The estimation of the equivalent height of the 
column involves a two-fold assumption ; first an 
estimation of the surface of the liquid on the top 
of the column and second an estimation of the 
velocity at the surface. Both are necessary in 
order to estimate the time of contact of an 
element of surface area. However, the value of 
Ah amounts to only 5-10%, of h, so that a relatively 
large error in the estimation of Ah will not greatly 
alter the relationship of the gas absorption to 
height, since this varies with h,}. 


The estimation of the actual temperature of 
the interface is subject to the most uncertainty 
of all the assumptions involved. If the film 
actually behaved as a semi-infinite stagnant 
layer, the surface temperature would remain 
constant. However, the thermal diffusivity of 
water is of the order of a hundred times greater 
than the diffusivity of SO, in water, so that, 
whereas the assumption of constant velocity and 
infinite depth is justifiable for the diffusion of 
material in the films considered, it is not at all 
justifiable for the diffusion of thermal energy. 

One can calculate the temperature rise of the 
surface of the film due to the absorption of SO, 
(Danckwerts [14]). If one takes the heat of 
solution of SO, in a saturated solution to be 
6-7 Keal/mol (Roru [12]}) this rise is found to be 
about 10°C. For want of a better assumption 
the average temperature of the surface of the 
film was taken as the average of the inlet and 
outlet water temperature plus 1°C. 

It is estimated that the total uncertainty due 
to systematic errors is + 5%, the greatest part 
being due to the temperature of the surface of the 
film. 


6. Limits oF APPLICABILITY 


The conditions studied here were for effective 
column heights between 12-5 and 22-5cm and for 
Reynolds’ numbers between 123 and 1150. This 
corresponds to contact times between 0-13 and 
0-98 seconds. Lower flow rates were not studied 
because of the difficulty in keeping the tube com- 
pletely wet by the solution. Measurements with 
shorter tube lengths were not included because of 
the increasing relative importance not only of the 
entrance effect but also of an exit effect. 

This exit effect is evidenced by the appearance 
of a band of ripples about a centimetre above the 
level of the liquid in the receptor. Preliminary 
qualitative experiments indicate that the surface 
of the falling film below this band of ripples is 
nearly stagnant, and that the height of the band 
of ripples varies inversely wih the velocity of the 
falling film. A similar phenomenon has been 
reported by Marsuyama [15], who studied 
absorption by short liquid jets. This effect will 
be discussed in Part II of this paper. 


The authors wish to acknowledge the contribu- 


Absorption studies in the light of the penetration theory - I 


ACKNOWLEDGMENT Ah = increment in height to account for entry effect, em 
h, = effective height of column, h + Ah, em 
r = radius of tube, cm 


tions to this paper by A. Kessier, who developed Re = Reynolds’ number of film 
the absorption apparatus, and to R. A. T. O. T = absolute temperature, “K 
NiJsinG, who helped with part of the taking of t = time, see 


the data. {, = contact time of an element of area in the column, 
sec 
vj = velocity of the surface of the liquid film, em/sec 
NotaTiox z = coordinate in the direction of diffusion, em 
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Absorption studies in the light of the penetration theory 
II. Absorption by short wetted-wall columns 


S. Lynx, J. R. Srraarewerer and H. Kramers 
Laboratorium voor Physische Technologie, Technische Hogeschool, Delft 


( Received 15 November, 1954) 


Summary~ The rate of absorption of SO, by water in short wetted-wall columns (heights 
of 1 to 5m) with no gas phase resistance has been investigated. When the liquid film emerges 
through a slit of appropriate width, the entrance effect is negligible and need not be considered 
in applying the penctration theory to the rate of absorption. However, there is a considerable 
end-effect which is evidenced by a thickening of the film, accompanied by an apparent stagnation 
of the surface, on the lower one centimetre of the column. The rate of absorption by this part 
of the film is relatively low. This phenomenon complicates the analysis of data obtained with 
short wetted-wall columns, or other equipment with films flowing over small distances. 


Résumé—On a mesuré la vitesse d'absorption de SO, pure dans eau dans un absorbeur a 
paroi moillée avec un hauteur de laSem. Si Peau est alimenté par une fis-ure d'une largeur 
appropriée, la correction d'entrée, qu'il faudrait introduire dans la theorie de la “ pénétration ™ 
pour calceuler Pubsorption, parait étre negligeable. autre part, il se produit une région avec 
une grande resistance pour labsorption a la partie inférieure de la surface de l'eau. En certains 
cas cette partic peut avoir une hauteur de 2cm. L‘apparition de cet effet est visible par une 
faible augmentation de lépaisseur de la couche, qui est accompagné par une corrugation de la 
surface. Par leffet signalé ici analyse de la résistance de la phase liquide & l'absorption se 
complique, surtout quand il s‘agit d'un appareil, o8 des couches de liquide parcourent des trajets 


courts. 
1. INTRODUCTION t, = time during which an element of area 
In the first article of this series [1], the absorption on the surface is exposed to gas (sec) 


of SO, by water using relatively long wetted-wall 
columns was discussed. Columns with effective 
heights between 12-5 and 225m were studied 
using pure SO, (thus no gas-phase resistance), and 
were found to obey the penetration theory quanti- 
tatively. It was shown that the film behaved as 
a semi-infinite layer with the surface suddenly 
saturated with SO,, and that the rate of absorption 
®,, (g/sec), for the column was expressed by the 


In the previous study the water flowed down 
the outside of a vertical tube which was open at 
the top. This produced an entry effect whose 
magnitude could not be estimated with great 
accuracy. For the study of short columns the top 
was covered with an accurately centred cap, which 
formed an entry-section for the film of water. The 
distance required for the surface of the film to 
; accelerate to 90°, of the equilibrium velocity was 
following formula : —_ estimated by the use ea sealed-up model and 

@, = 2(C* —C,) A af 2. (1) found to be about 12 times the film thickness. 
#G Thus, the entry-effect is negligible and need not 
where C* = concentration of solute at the sur- be considered in applying the penetration theory 


face (g/cm*) to the rate of absorption. 
C, = initial concentration of solute in bulk In the course of the experiments with both the 
of liquid (g/cm*) long and the short columns, it was found that at a 
A = area of film surface (em*) certain height above the surface of the stagnant 
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liquid in the receiver, into which the film of water 
flowed, the thickness of the film increased from 
the normal to a slightly greater value. This was 
accompanied by the appearance of a narrow, 
horizontal band of ripples at the top of the 
thickened portion of the film. The height of this 
end-effect was found to vary in an inverse manner 
with the velocity of the film, from 2m for the 
lowest flow rate used to 0-1 em for the highest. An 
effect similar to this has been observed by 
Matsuyama [2], who studied absorption by a 
short jet of liquid. However, the effect apparently 
escaped the attention of Praceman [3] and 
Ewvmenrt and Picrorp [4], who have also worked 
with short columns. 

Qualitative experiments indicate that the 
surface of the film below the ripples is nearly 
stagnant. A light powder (lycopodium seed) blown 
onto the surface of the film circles lazily instead 
of being swept immediately into the liquid 
receiver. A coloured liquid introduced in the 
surface layers of the film accumulates to form a 
dark-coloured band below the ripples [5]. 

If the surface of the film below the ripples were 
actually stagnant, the rate at which that area 
absorbed gas would be very much lower than that 
of the much fresher liquid above the ripples. One 
may consider two extremes: (1) when the area 
of this end-effect is a comparatively small fraction 
of the total film area and (2) when the end-effect 
has reached the top of the column so that the 
surface of the entire film is apparently stagnant. 

In the first situation, the rate of absorption 
of the area below the ripples may be neglected 
relative to that of the fresher surface. The column 
may thus be considered to be shorter by the height 
of the end-effect, h, Using the relationship 
between the contact time, ¢,, and the liquid flow 
rate, acceleration of gravity, and the physical 
properties of the fluid given in [1], equation (1) 
then becomes 

3 
= 4(C* —Cy)(r + Ph, (2) 
where h, = effective height of column, (h — h,), 
cm 


In the second situation it is assumed that the 
liquid issuing from the entry-section, remains in 


full parabolic flow, with zero velocity at the 
surface. It is thus analogous to the case of heat 
transfer to a liquid in laminar flow between 
parallel plates when the liquid enters a section 
where one plate has a different, constant tempera- 
ture and the other is perfectly insulated. So far 
as is known, the literature contains no solution 
of the differential equation obtaining in this case. 
However, for sufficiently low values of the term 
hD 5T., the solution may be approximated from 
the solution for the symmetrical case when the 
liquid enters a section where both plates have the 
same different, constant temperature. The rate 
of heat transfer in the former case will be just 
half that in the latter case, in the first approxima- 
tion (in reality it will be somewhat greater than 
half). The solution for the latter case has been 
calculated by Prins, Mutper and Scuenk [6]. 
Using the notation of this paper, their results 
may be expressed in the following formula : 


where 5, for flow due to gravity between vertical 
parallel plates, is 


The value of the term 


(3) 


function of pa , calculated from the work of 


Prins et al., is given in the following table. 


Table 1. 
hD 
sl, 
0-01 0-025 
0-02 0-040 
0-05 0-075 
0-10 0-120 


2. Description OF EQUIPMENT 
The equipment has been described in detail in 
[1]. The chief modification for the purposes of 
this paper was the introduction of the water 
through a narrow slit at the top of the column. A 
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second alteration was the addition of a plastic 
cover to the liquid receiver. The liquid in the 
receiver was kept in. contact with the cover during 
a run so that the surface exposed to the SO, was 
greatly reduced. A sketch of the absorber, showing 
the changes mentioned, is given in Fig. 1. For 
further details of the equipment the reader is 
referred to [1). 


WATER OUT twaree 


Fic. 1. Sketch of modified absorber. 


38. TREATMENT OF Data 


The data were reduced to the common basis of 
a partial pressure of SO, of 750mm Hg and a 
surface temperature of the water film of 20°C. 
The manner in which these adjustments were 
made has been described in [1]. In all experiments 
pure (99-9%) SO, was absorbed in distilled 
water. In some cases a surface-active agent 
(Teepol) was added. The height of the band of 
ripples was measured with a millimetre scale 
mounted inside the absorber with a precision of 
about half a millimetre. This measurement was 
made while gas was being absorbed. 


4. Discussion or REsvutts 


Gas absorption measurements were made at 
column heights of 5, 3, 2, and lem using pure 
water and dilute solutions of Teepol of various 
concentrations. These data are plotted in Fig. 2 
as ®, versus I+. For comparison, equation (2), 
with h, = 0, has been drawn in for the four 
column heights, and equation (3) has been drawn 
in for a height of 2cm, The values of C* and D 


used were taken from [1]. For the shorter 
columns equation (2) was drawn in only for values 
3 
of above 1-2 
em see. 
with regard to which points should approach 
which lines. 


to avoid confusion 


@ 4=30cm 02000 = 
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Fic. 2. Rates of absorption for different heights of column 
as functions of flow rate. 


In Fig. 2 the effect of the stagnation of the 
surface of the film below the band of ripples is 
quite apparent. As the flow rate increases, the 
height of the band of ripples decreases, and the 
measured rate of absorption approaches the 
calculated rate. 

The height of the end-effect is not affected by 
the height of the column (see Fig. 3). For this 
reason the fraction of the total area of the column 
which is below the band of ripples at a given 
liquid flow rate will be greater for the shorter 
columns. This is reflected in the rate of absorp- 
tion, as shown in Fig. 2. The deviation of the 
points from the theoretical lines shows a con- 
sistent increase as the column height decreases. 

Various Teepol concentrations were used with 
the 3-cm column. It is seen that, within the 
experimental error, there is no observable effect 
of the addition of the Teepol. 

For both the 1- and the 2-cm columns the band 
of ripples reaches the top of the column at a 
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certain flow rate below which the entire surface 
of the film appears to be nearly stagnant. The 
absorption data for the two columns show distinct 
breaks at those flow rates. 

Absorption data for the 2-cm column were 
taken at exceptionally low flow rates, in the 
region where the entire column is covered with the 
stagnant film. As seen in Fig. 2, the points lie 
consistently 15 to 20% above the curve of 
equation (3). It is interesting, however, that the 
agreement is as good as this considering the 
assumptions involved in the derivation of the 
equations. 

In Fig. 3 is a plot of +/h, versus I}. It is seen 
that the height of the end-effect does not depend 
upon the height of the column, but is influenced 
by the Tepol concentration. As Teepol is added, 
the band of ripples becomes fainter and more 
difficult to observe. At low Teepol concentrations 
h, is somewhat greater than with water, while at 
high Teepol concentrations it is appreciably 
lower. As mentioned before, however, the addi- 
tion of the Teepol had no measurable effect on the 
gas absorption. For this reason, the line shown 
in Fig. 3 was drawn through the points for water 
without regard to the points for the Teepol 
solutions. 

The straight line in Fig. 3 was used to correlate 
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Fig. 8. Height of the end-effect as a function of flow rate. 
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the data for various columns in Fig. 4, where ®,, 
is plotted versus I} (h — h,)!. Equation (2) 
is again drawn in the figure, and it is seen that 
the correlation is remarkably good. From the 
agreement of the data with the calculated line 


dh 
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Fic. 4. Correlation of rates of absorption for different 
columns by subtracting stagnant area. 


it must be concluded that the stagnant lower 
part of the film indeed absorbs gas at a greatly 
reduced rate. 

In Fig. 5 the data for the 1-cm and the 5-cm 
columns are plotted on a log-log scale to show 
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Fic. 5. Log-log plot of absorption data to demonstrate 
weakness of this form of presentation, 
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the interpretation that one might be led to from 
such a presentation. The points for the 5-cm 
column lie on a rather nice straight line having a 
slope of 0-39. If one discards two of the points 
for the l-cm column, the other six form a line 
witha slope of 0-97. The data for the other columns 
ean be plotted similarly to obtain values for the 
exponent n in the equation 


(4) 


which lie between the two mentioned above. It 
is possible that the presence of stagnant regions 
on the liquid films flowing in packed columns is 
responsible for the well-known facts that the 
exponent in equation (4) is usually about twice 
the value of 1/3 predicted by the penctration 
theory, and that the value of this exponent de- 
pends upon packing shape and size and the 
nature of the absorbing liquid. 

The presence of such stagnant regions can be 
demonstrated by introducing lycopodium powder 
into a laboratory-size packed column. Whirling 
groups of powder particles can then be seen, and 
the residence time of the powder in the column 
is found to be far greater than that of the liquid. 


5. Discussion or 
The work presented in this article is subject to 
the same errors and uncertainties discussed in 
[1]. Because the rate of gas absorption by the 
short columns is smaller than by the long columns, 


the percentage deviation of the experimental 
points from a smooth curve drawn through them 
is somewhat greater in this work. 


6. Limits or APPLICABILITY 


The conditions studied here were for column 
heights between 1 and 5cem and for Reynolds’ 
numbers between 20 and 1150. This corresponds 
to contact times between 0-01 and 0-22 seconds 
neglecting the stagnation effect. 
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NOTATION 


A = urea of film surface, em? 
Co = initial concentration of solute in bulk of liquid, 
g/em® 
C* = concentration of solute in saturated solution, g/em*® 
D ~ diffusivity of solute in liquid, em? /sec 
g = acceleration due to gravity, 981 cm /sec* 
h = height of column, em 
h, = height of band of ripples, em 
h, = effective height of column (h — h,), em 
r = radius of tube, em 
Re = Reynolds’ number of film 
{. = contact time of an element of area in the otumn, sec 
I’, = volumetric rate of flow per centimetre, —-— 
em /sec 
8 = thickness of liquid film on wall of column, em 
vy = kinematic viscosity, em? /see 
®,, — rate_of absorption, g/sec 
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Absorption studies in the light of the penetration theory 
III. Absorption by wetted spheres, singly and in columns 
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Summary—A study has been made of the absorption of SO, by water flowing in a laminar 
film over the surface of single spheres and columns of spheres. By means of the penetration 
theory it is possible to correlate the data for spheres of radii from 1 to 3. cm, and for columns 
containing up to fourteen spheres. It appears that there is little or no mixing of the surface 
layers with the bulk of the liquid as it flows from one sphere to the next in a vertical row. 


Résumé—On a é¢tudié l'absorption de SO, par eau, qui coule en couche laminaire sur la surface 
d'une sphére et d'une colonne formée par sphéres superposées. I] a paru possible de trouver 
une correlation satisfaisante entre la vitesse d’absorption et la dimension des sphéres (2 a 6 cm. 
diamétre), leur nombre (1 4 14) et le debit d'eau (0-5 a 7-5 em*/sec). Cette correlation peut étre 
prédite au moyen de la theorie de diffusion non-stationnaire. Les expériences démontrent qu’il 
ne se produit pas un melange entre la surface et l'intérieur du liquide quand celui-ci passe d’une 


sphére a l'autre. 


1. INTRODUCTION over the sphere in a laminar film, with a total 


In the first article of this series [1] it was found Volumetric flow rate of ®,(em*/sec). It is 
that for physical absorption without gas phase assumed that the thickness of the film, 5, at any 
resistance in a wetted-wall column the simple latitude on the sphere is the same as it would be 
theory of transient diffusion into a stagnant layer for the same flow rate per unit length, I (em®/- 
of liquid (penetration theory) is valid for contact ©™-sec), on a plane surface making the same angle 
times up to 1 second. In the case of a film of With the vertical. Thus : 


liquid flowing over a sphere a similar mathematical 4/30. 

treatment is greatly complicated if one attempts 8 = J sane 

to consider the stretching of the film, and the oT 

consequent distortion of the concentration = Sint (1) 


gradients, as it flows from the top of the sphere 
to the bottom. This stretching will have the effect where » = kinematic viscosity, em?/sec 

of increasing the concentration gradients on the g = acceleration due to gravity, 981 em/ 
upper half of the sphere and of decreasing them nec®, 

on the lower half. Because the film on the upper 
half of the sphere, being fresher, absorbs faster 
than when it is on the lower half, one would 
expect the net effect of the stretching and con- 
tracting of the film to be an increased rate of 
absorption. Without considering these effects, 


It follows from the above assumption that a half- 
parabolic velocity profile will exist at all latitudes 
on the sphere, so that the velocity of the surface 
of the film, v; (em/sec), is given by the following 
equation : 


the following calculation can be made : » = 8 

Consider a sphere of radius R(cm) with an 
angular coordinate « which is zero at the top of og \* © (2) 
the sphere and at the bottom. Liquid flows (ax: a) 
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The rate of absorption of a differential unit of 
area on the film is 


d®, =(C* —C,) [P aa, 


where @,, = rate of gas absorption, g/sec 
C* = concentration of solute in saturated 
solution, g/em® 
C, = initial concentration of solute in bulk 
of liquid, g/em® 
D = diffusivity of solute in liquid, em?*/- 
sec 
t = time during which the element of 
area has been exposed to the gas, sec 


The “age”’ of an element of area at latitude « is 


(8) 


(5.7) ‘0.2/9 sin’ ada 


and an element of area can be expressed as 


dA = R*sinada«. 


Thus, the rate of absorption for the whole sphere 
becomes 


©, = (C* —C,) VD 


sinada 


fai 


(4) 


The value of the term 
"sinads i 
in’? ada 
and the value of 
fine x is 2-58, as 


determined by numerical integration. 


For the absorption of SO, in water when the 
partial pressure of the SO, is 750mm Hg and 
the temperature of the surface of the water film 
is 20°C, equation (4) can be expressed as 


,, = K (4a) 


The value of K obtained using the values of 
the pertinent physical properties listed in [1] is 
1-19 x 10% [g SO, /sec (em* H,O/sec)”*). 
This value can only be considered as an approxi- 
mation because the flow conditions at the top 
and bottom of the sphere have been idealized and 
the effect of the stretching and contracting of the 
film has been disregarded. 

The derivation above has been for single spheres. 
For spheres mounted in a column one must con- 
sider whether the streamlines of the flowing film 
are effectively mixed as the liquid passes from 
one sphere to the next. As has been shown above, 
the rate of absorption in the case of one sphere is 


A 
~(C* — O) i, 


where ¢, is the time in which an element of area 
is exposed to the gas in moving from the top of 
the sphere to the bottom. If there is complete 
mixing of the streamlines of the fluid at the 
points between the spheres, and if the concen- 
tration of the dissolved gas in the outgoing liquid 
is small relative to C*, the rate of absorption for 
N spheres will be approximately : 


] spheres =N- [,,]1 ephere 
If there is no mixing of the streamlines in the 
flow from one sphere to the next, the contact 
time for the liquid flowing down a row of N spheres 
will be N -t,, the total area exposed to the gas 
will be N - A, and the rate of gas absorption will 


be 
or from equation (4a), 
= K /N (5) 
It should be remembered that the derivation 
of these equations rests on the assumption that 


film can be considered as a semi-infinite medium. 
With other words, the “ depth of penetration ” 
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of the solute must be small, and the relative 
saturation of the outgoing liquid will be low. 


2. Description or EQUIPMENT 


The equipment used was the same as that des- 
cribed in [1] and [2] except that the absorber 
was modified to permit the placement of the 
spheres, as shown in Fig. 1. Single spheres having 
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Fic. 1. Sketch of modified absorber. 


radii of 1-00, 1-50, 1-90, and 2-95 cm were used. 
Vertical rows of 3, 5, 10, and 14 spheres having 
radii of 1-00 cm were also studied. The bottom 
of the lowest sphere was set 2 cm above the level 
of the liquid in the receiver to minimize the in- 
fluence of the end-effect described in [2]. The 
liquid receiver was equipped with a plastic cover, 
also described in [2]. For further details of the 
equipment the reader is referred to [1]. 


8. TREATMENT OF Data 


The data were reduced to the common basis of a 
partial pressure of SO, of 750 mm Hg and a sur- 
face temperature of the water film of 20°C, The 
manner in which these adjustments were made 
has been described in [1]. In all experiments pure 
(99-9%) SO, was absorbed in distilled water to 
which a surface-active material (Teepol) had been 
added. The spheres were made from steel ball- 
bearings and were mounted on a stainless steel 
rod. 


The results presented here are subject to the 
same errors and uncertainties discussed in [1] and 
[2]. The average deviation of the points in Fig. 2 
from the solid line is 0-11 10% g SO, /sec. 


4. Discussion oF RESULTS 
The rate of absorption using single wetted spheres 
is shown in Fig. 2, plotted against +/N R7/* 61/8. 
The slope of the solid line through the points is 
about 7-6°%, greater than the value of K calculated 
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Fic. 2. Gas absorption by single spheres. 


from the penetration theory. This disagreement is 
probably due to the stretching effect, which was 
neglected in the derivation of equation (5), and 
is actually smaller than might have been expected. 
Of more importance is that the data for the four 
spheres all fall on the same straight line when 
plotted in the manner suggested by the pene- 
tration theory. 

The solid line of Fig. 2 is also drawn in Fig. 3, 
where it is seen to represent the data for vertical 
rows containing up to 14 spheres with good 
precision at sufficiently low flow rates. It is seen 
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Fic. 3. Gas absorption by columns of spheres, R = 1-00 em 


that the data for the 5-sphere column begin to 
deviate appreciably from the solid line at a flow 
rate of 4:2cm’/sec. Those for the 10-sphere 
column begin to deviate at 2-0 cm*/sec, and those 
for the 14-sphere column begin to deviate at 
1-9 cm®/sec. This effect is due to the increasingly 
great tendency of the liquid film to ripple at a 
given liquid flow rate as the number of spheres 
in the column is increased. The instability of the 
flow over the lower spheres in a column is ob- 
served as a tendency for the liquid stream to 
“dance” back and forth at the points where 
the spheres are joined to each other. This tendency 
is decreased, but not eliminated, by the addition 
of 0-050%, of a surface-active material (Teepol). 
With concentrations of Teepol either greater or 
less than this the tendency to ripple is greater. 

The fact that the data for the columns of 
spheres fall on a straight line when plotted in 
the manner of Fig. 3 is a strong indication that 
there is essentially no mixing of the streamlines 
in the film as it flows from one sphere to the 
next. 


These results suggest that liquid in laminar 
flow in a packed column may frequently travel 
over many pieces of packing before effective 
mixing of the surface layers and the bulk of the 
liquid occurs. The average liquid-side resistance 
to absorption would then be greater than if 
mixing were taking place at each juncture be- 
tween the pieces of packing. However, as has 
been indicated in [2], at the places where such 
mixing does occur, there is a stagnant surface 
over the region one to two cm above the point 
of mixing, again resulting in a reduction of the 
absorption rate. In the practical case, the effects 
described in this paper and in [2] are counteracted 
by turbulent regions in the film, which are pro- 
duced by local uneven wetting of the packing 
and by the stream of gas. This turbulence tends 
to decrease the liquid-side resistance. 


5. Limits oF APPLICABILITY 

The conditions studied here were for single 
spheres with radii between 1-00 and 2-95 em, for 
columns containing up to 14 spheres having radii 
of 1-00 cm and for flow rates between 0-5 and 
75cm? H,O/sec. This corresponds to contact 
times between 0-29 and 0-05 sec for the 1 cm- 
radius sphere and between 1-78 and 0-29 sec for 
the 3 cm-radius sphere. 


ACKNOWLEDGMENT 
The authors wish to acknowledge the contribu- 
tions to this paper made by A. J. pe Graar and 
J. C. Hoocrveen, who took part in the investi- 
gation. 


NOTATION 
A = area of sphere, em? 
C* = concentration of solute at saturation, g/em*® 
C, = initial concentration of solute in bulk of liquid, 
g/em* 
D = diffusivity of solute in liquid, em? /sec 
g = acceleration due to gravity, 981 em /sec* 
K = constant in equation (4a) 
N = number of spheres in a column 
R = radius of sphere, em 
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t = time of contact with gas of an element of area, sec ‘J’, = liquid flow rate per unit length, cm*/cm sec 
t = time of contact with gas of an element of area in 5 = film thickness, em 


moving over one sphere, sec v = kinematic viscosity, em? /sec 
0; = velocity of surface of film, cm/sec ®,, = rate of absorption of gas, g/sec 
a = angular coordinate on sphere ®, = liquid flow rate, em*/sec 
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On the steady state fractionation of multicomponent and complex 
mixtures in an ideal cascade 


Part 2.—The calculation of the minimum reflux ratio 


Axpreas Acrivos* and Neat R. Amunpson** 
Department of Chemical Engineering, University of Minnesota, Minnneapolis 14, Minnesota 


( Received 6 October 1954) 


Summary = The important problem of calculating the minimum reflux ratio in the separation 
of multicomponent and complex mixtures in an ideal cascade has been solved in closed form. 
The solution of two rather unusual integral equations may be of interest. Numerical calculations 
have been performed and will be reported elsewhere. 


Resumé — Les auteurs ont résolu le probleme important du calcul du rapport du reflux minimum 
dans la séparation d'un mélange de composants multiples et complexes, dans le cas d'une cascade 


idéale 


La solution de deux équations integrales assez particuliéres peut étre intéressante. 


Les calculs numériques ont été faits et seront présentés dans un autre texte. 


5. INrROopUCTION 


In a previous article the authors examined 
theoretically the fractionation of multicomponent 
and complex mixtures in an ideal cascade and 
solved the equations, describing the process, in 
a closed form. The object of such a calculation 
is the following. Given a mixture of known com- 
position and thermal condition it is desired to 
determine the number of stages in the enriching 
and stripping sections of the cascade to separate 
the feed into two streams such that 2,(8) and 
are specified a priori. Here 8 and are 
two components and «,(j) is the mole fraction of 
the jth component in the product streams. How- 
ever, the final answer depends on two additional 
parameters, the feed stage location and the reflux 
ratio R. One introduces the feed at such a 
location to minimize the total number of stages 
required for the specified separation. It is found 
that as the reflux ratio is decreased the number 
of plates in both sections of the column, N and 
S, increase until a value of R is reached for 
which both N and S are essentially infinite, even 
with the optimum feed stage location. Since in 


an actual operation R > R,,, the calculation 
of this minimum reflux ratio for any desired 
separation is of considerable importance. 

Exact methods for determining R,,,, in an 
ideal cascade have been worked out recently for 
discrete mixtures by UNpERwoop [1] and 
Murpocn and [2]. In this paper the 
formulae of Murxpocu and will be 
derived in a somewhat more general form, and 
these will be generalized to the complex mixtures 
when the number of components in the mixture 
becomes infinite. 


6. Minimum Reruux ror Discrete 


MIXTURES 

Equation (3.19) with p,* = p,; may be written 
a,(i) R 

c,* 

pt) (6.1) 

)N-1 
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Let the two functions 
hn (6.2) 


be defined where | and fA are arbitrary for the 
moment save that 1 < h, 1 < m. Then (6.1) can 
be written 


m G,* 


R* P(t) — doe (6.4) 
ro*(i) R > G, 
kat Ay — 
and 
Z 
G,° 


Now the numbers 4,, 1 < k S m, are the m roots 
of the characteristic equation (8.7), i.e., 


PJ) _ 

and therefore one can say 
P(A) = 


A- PAI) Fo\J) = 


From this it follows that if A = p(i) then* 


— pli) 
ont (6.6a) 
(p(k) — pli)) 
Similarly 
(Plt) — $e) (6.6b) 


(pli) — plk)) 
kel 


Finally if r;,i sj S m—1, are the m — 1 finite 
roots of 
Ae 7; 
* a, indicates an obvious term must be on.itted from the 
product. 


=0 (6.7a) 


it can be shown that 


k=l 

(Ae — P) 
and 


w (p—r;*) 

-= See (sp) 


where r;*, 1 Sj S m—1, are the m— 1 finite 
roots of 
m G,* 
= 
— 
Therefore if the above equations are substi- 
tuted into (6.4) it is found that 


7 (r,* —p,),1 


0 (6.7b) 


and hence 

r, 1sksm-1l1 (6.9) 
The function G, will now be examined for large 
N. By definition 


p(m) < p(m — 1) < ---< pili) 


< A < p(k — 1) 
Hence if it is assumed that C, + 0, then since 


N-1 
it follows that as N + 
(a) k>h 
(b) G, +> 1 by definition 
(c) for k< h only it tim 


But this implies that fork < h,C, +0asN + @ 
and therefore 
P(k) 
* 
Hence it is included that G,, k < h, exists, if and 
only if 
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tim (6.10) 
An 
exists. Moreover, it is seen that as N + o@ the 
function 1,(i) can be separated into two distinct 
parts such that 

(a) a,(i) 2h 
(b) agli) +0, i < h, as N + @ in such a way 
that 
A.\* 
exists. Similarly 
(a) +0, i sl 


(b) +0, i > 1, as S + in such a way 


that 
lim a,*(7) 
¢, 
exists and G,* +0 fork = l, while G,* + 0, for 
k 2 l as S +o. 

It follows then from (6.7a), (6.7b), and (6.9) 
that as N + wo, S > o 

and, therefore, at minimum reflux 


The above result was first derived rigorously by 
Murpocu and Hoiianp [2] although it had been 
anticipated earlier by UNpexwoop [1]. 

However, the two characteristic equations from 
which A, and ¢, are calculated are 


and since 

+ Wag(j) = 


where 2,(j) represents the mole fraction of the 
jth component in the feed, it follows from the 
above and from (6.11) that 


Pj) 
— PG) 
The parameter g = L* —L can be calculated 


(6.11) 


(6.12) 


—q-h+1 sk sl (6.13) 


from a knowledge of the thermal condition of 
the feed. 

It is seen then that (1 — A) characteristic roots 
Ap (kh +1) 5 SI can be determined from 
(6.13) without a knowledge of the function z,(j). 
Also, since x,*(j) = 0 for j > I, it follows from 
(6.11) that 

j>l (6.14) 
Now the characteristic equation for the enriching 
section of the cascade is 

i . . 

Pj) Pj) Daj) 

since 2,(j) = 0 for j < h. 

Next two functions 4, and B(k) are defined 

such that 
A 
= 
Ay — Pli) 
and 
Blk) 
Aj Pk) 
It can be shown easily that 

h+1sijsl 


1, h+1sisl (6.16) 


h+1<i sl (6.17) 


=1, 


(6.18) 
and, moreover, that 


(Ay — pli)) 
Bi) = -,h+1sisl (6.19) 
(p(k) —p(i)) 
So by multiplying (6.15) through by 4,/(A,—p(i)) 
and summing on k there results, because of (6.18) 
and (6.14) 


p(h) D xo(h) 
Pi) = Bl) [RD ~~ Bh) (pli) — p(h)) 


forh + i Sl where 


— 


(p(k) — p(i)) 
i=handizl+1 


| (6.20) 


(6.19) 
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In this manner (/ — h) linearly independent non- 
homogeneous linear simultaneous equations in 
(1 — h) unknowns have been solved. It must be 
kept in mind that A,,h + 1 Sk <1, and B(i) can 
be calculated a priori from the known composition 
and thermal condition of the feed. Hence from 
(6.20), from 


Dz{j) = (6.21) 


j=<h 

and the fact that 2,(8) and 2,(8*), where h = 8, 
8* <1 are specified, one can calculate both the 
complete composition of the overhead product 
and the minimum reflux ratio. Of course, it 
must be understood that A and | are not known 
a priori, and so, in general, there are more 
unknowns than there are equations. However, it 
turns out, in all the cases which have been investi- 
gated so far, that the system of equations (6.20) 
and (6.21) has only one solution which satisfies 
the additional requirement 


0< Da,fi)<2,fi) for (6.22) 


In practice it is best to assume A and / first, and 
then to solve equation (6.20). This process is 
repeated until (6.22) is obeyed. 

The solution to the minimum reflux problem 
is exact. It is desired now to develop the analogous 
results for complex mixtures. 


7. Minimum ReEFiux For Complex 
MIXTURES 


1. Continuous mixtures. By definition the dis- 
tribution function of a continuous mixture is 
single valued and continuous and it will be 
supposed that 2,(p), the distribution function of 
the overhead product, does not vanish along the 
positive p axis. Therefore from (4.26) 


ay(p) = 
N-1 C(p) C(p) 
(7.1) 


Or, if the new function G(A) is defined 
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C(A) 
00) = (i ) 
where A is, thus far, arbitrary, then 
= 
G(p) 
Pi —— da + 
A—p PtP) 
(7.8) 
| GA) da 
Now as N > o, it follows from (3.2) that if 
C, +0 
(a) GA)>0, A<h 
(b) G(h)=1 
C(A) 
ly if | 
(c) > h, exists only i im (; 


But, according to the definition of C(A) 
A A) 
CO) = + 


and therefore it may be concluded that G(A) for 
A > h exists only if 


A) 
lim (;) 
Now h 
exists. 
Some complications arise, however, because as 


N becomes large and is made formally to approach 
infinity, one obtains from (7.3) 


GA) 
pf 


should exist which implies that 2,(h) and therefore 
C(h) approach zero as N + o. In spite of this, 
it can be rigorously shown that 

lim G(A)=0 , A<A. 

Thus, at minimum reflux the distribution function 
of the overhead product, 2,(p), can be separated 
into two distinct parts, such that 


(a) a(p) +0, p<h 
(b) a(p) > 0, p > h, according to (7.5) 


(7.4) 


(7.5) 


lim a,(h) (7.5a) 
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while z,(h) approaches zero according to (7.5a). 
Similarly the distribution function z,*(p) can be 
separated into two distinct parts, such that 


(a) + p>l 
(b) 2,*(p) > 0, p < J, in such a way that 


To solve the problem for continuous mixtures, 
one sub-divides the interval 0 < p < a, where a 
is arbitrary, but large, into m equal parts, Ap, 
and determines the asymptotic form of the 
solution for the discrete problem as Ap + 0. The 
limiting form of (6.13) is 


(7.5b) 


for | = A < h, an equation which defines a func- 
tion p(A), | SA< A, which is identical for, 
1s A< A, with the function p(A) defined by 


| ap + Aag(A) weot mp(A) = (7.6a) 


On the other hand (6.19) becomes 
p(k) 


= ——-,, OF 


Ap 


But since, 0 < p(k)< 1, the above can be 
expanded in a Taylor series, and it is found that, 
asm, > @ 


B(p) = 
p(p) exp J da — (7.7) 
1 


In = In pli) + In 
kel 


for 1 <= p< h, and where Z(2i) is the Riemann 
zeta function. Similarly it can be shown that 
(6.19a) becomes 


B(p) = exp (fea) aa) 


for p < 1. It is worth pointing out that (7.7) is 
the solution to the integral equation 


B(p) 


ls A< h, where p(A), 0 < p(A) < 1, is a known 
function. It can be shown that as p > /| 


Bip) + K 


and therefore B(p) is integrable. 
Finally, the limiting form of (6.20) becomes 


i 
q2,(q) 


D = RD 
p Dap) + PONG 


for! < h, while 
RD 
B(q) — 


since B(l) is infinite, and 


D 


(7.9a) 


(7.9b) 


It is again interesting to point out that (7.9) is 
the solution to the integral equation 


+ A Da,(A) cot  p(A) 


forl SA<h. 

Thus far it has been assumed that z,(p) did 
not vanish along the p axis. Let it be supposed 
now that z,(p), | <p Sh, vanishes only and 
everywhere in the set of intervals B, b, < p< a, 
k=1,2,...¢, where a, >a,... >a, and let 
the remaining part of the interval 1 <p <h 
be A. 

Now the limiting form of (6.18) is given by 
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(7.6) for A in A, but, in addition, c numbers 
Ay A, in B, can be found such that 


(7.10) 


It follows then that the limiting forms of (6.19) 
and (6.19a) are respectively 


— 
= exp 
(by — P) 
n 
1 


for all p in A, and 


— p) 
Bp) = + 
(b, — p) 


(7.11a) 


for p < l. With these new definitions Eqs. (7.9), 
(7.9a) and (7.9b) remain unchanged. 

It must be kept in mind, however, that the 
two points h and / are, in general, not specified 
a priori. But, at any rate, it is best to assume h 
and / first and then to solve (7.9) and (7.9a). This 
process should be repeated, until the calculated 
values for 2(p) at the two points p= 8 and 
p = B* coincide with those stated in the problem. 

2. Semi-continuous miztures. By definition the 
distribution function of a semi-continuous mixture 
is everywhere single valued and continuous except 
at a finite number of points where it possesses 
infinite but integrable discontinuities. In section 4 
of the previous paper, the fractionation of semi- 
continuous mixtures in an ideal cascade were 
considered. The minimum reflux problem for 
semi-continuous mixtures will not be worked out 
here except to say that the problem can be solved 
without much difficulty in the manner shown in 
the previous part of the present paper. It should 
be pointed out, however, that for every integrable 
infinite discontinuity in z,(p), for! << p < h, there 
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is an integrable infinite discontinuity in B(p) - 
defined by a proper generalization of (7-11) — at 
precisely the same values of p. 


8. CONCLUSIONS 


The minimum reflux problem in the rectification 
of multi-component and complex mixtures in an 
ideal cascade has been solved exactly. By means 
then of the formulae presented in this paper it 
is possible to calculate, not only the -minimum 
reflux, for a specified separation, but also the 
composition of the product streams. 

One of the main characteristics of this problem 
is the appearance of the two points h and | which 
are unambiguously defined. The value of both h 
and / will depend, of course, not only on the 
composition and the thermal condition of the 
feed, but also on the desired separation, which is 
specified by giving the value of 2,(p) at two 
points, say p= 8 and p = £*. In any actual 
calculations it is however preferable to assume 
h and I, until the calculated values of 8 and £* 
coincide with those given. 

In addition, the form of (7.3) and the analogous 
expression for the stripping column correctly 
predict that, in general, the cascade can be 
broken up, at minimum reflux, into five sections. 
The two sections, the so-called “ pinches,” which 
contain the infinite number of plates are, respec- 
tively, above and below the section containing 
the feed plate. Moreover, the composition of the 
various streams does not change, from plate to 
plate, in the two pinched reagions, whereas it 
does change in the other three sections. It can 
easily be shown from (7.3) that, in the pinched 
region of the enriching section, 


= 
for p < h, while x,(p) = 0 for p > h. Similarly, 
in the pinched region of the stripping section, 


— P%*(P) 


for p > l, while 2,*(p) = 0 for p < 1. The above 
relations hold for all mixtures. 


dp=-—q 
= 


Anpreas Acrivos and Neat R. AMUNDSON 


NoTaATION 
A, = function defined by Eq. (6.16) 
Bik) = function defined by Eq. (6.17) 
Cy = function defined by Eq. (3.11) 
C,* = function corresponding to C, for stripping section 
G, = function defined by Eq. (6.2) 
G,* = function defined by Eq. (6.3) 
h, | = key components 
N = number of plates in rectification section 
p(i) = reciprocal relative volatility of ith component 
q = thermal factor for column feed 
R = rectification section reflux ratio 
R* = stripping section reflux ratio 
1; = function defined by Eq. (6.7a) 
7;* = function defined by Eq. (6.7b) 
S = number of plates in stripping section 
Zo(i) = overhead composition function 
Zo*(i) = residue composition function 
Z = function defined by Eq. (6.4) 
¢, = characteristic root for stripping section 
Ay = characteristic root for rectification section 
p(A) = function defined by Eq. (4.15) 1955 
{(z) = Riemann zeta function 


Other symbols defined in Part I. 
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Thermodynamic evaluation of binary vapour-liquid equilibria 


H. A. C. Trisssen 
Laboratory of Physics and Meteorology, Agricultural University, Wageningen, The Netherlands. 
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Summary—An exact graphical method for the evaluation of vapour-liquid equilibrium data in 
a binary mixture from boiling point composition data is described. The variation of the activity 
coefficients with composition is obtained by successive application of the exact Duhem equation 
to a graph relating In y,, and In yg with 2 as parameter. Further, the test for internal consistency 
of vapour-liquid equilibrium data as introduced by Harincron is extended to p constant, T 
variable and corrected for nonideality of the vapour phase. The methods are applied to vapour- 
liquid equilibrium data of the binary system methylcyclohexane-toluene. 


Résumé- -L‘auteur décrit une méthode graphique exacte pour vérifier & partir de la courbe 


débullition d'un mélange binaire, les données de l'équilibre liquide-vapeur. Il obtient la variation 
des coefficients d’activité avec la composition par application de la forme exacte de l’équation de 
55 Dvunem a un graphique reliant In et In avec x comme parameétre. 


En outres, il étend le test introduit par Hartncron qui apprécie la qualité intrinséque des 
données de l’équilibre liquide vapeur, & p constant et T variable et le corrige dans le cas de non 
idéalité de la phase vapeur. I] applique cette méthode aux données de l'équilibre liquide-vapeur 
du systéme binaire methyl-cyclohexane — toluéne. 


Tue Gisss-DunemM Equation valid as equilibrium relation, since according to 


The general form of the Gibbs-Duhem equation the phase rule it is impossible to maintain 

for binary mixtures valid at variable T and p is ©4Uilibrium if the composition is varied at con- 
stant temperature and pressure. 

nm, dp, + ng du, = Vdp — Sat (1) F. DonGe [2] showed by a comparative 


If we replace the chemical potentials yz by the study of the convential Gibbs-Duhem equation 


fugacity ling to the definition of the and his derivation of the exact relationship valid 


fugacity at p constant and T variable that even for rather 
close-boiling mixtures the degree of inaccuracy 
du, = RTd\nf, at T constant in the extremes of the concentration range may 


t to 10° 
and divide by the sum of the mole numbers, 
eq. (1) te fom forins of eq. (83) as developed by van Laar, 
z,dinf, + = Oat Tandpconstant (2) Marcutes, Scatcnarp and others is that only 


two experimental constants, e.g. A, and A, are 
Introducing the partial derivatives of Inf with i, volved [1, 5, 9]. From a number of experi- 


respect to the mole fraction of one of the com- jnental vapour-liquid data an average A, and Ay 
ponents, eq. (2) reduces to the common Duhem oan thus be determined, and assuming eq. (3) to 


equation : be valid, an accurate calculation of the equili- 
brium curve can be made. 
= C24) or, since = fi 
T, Pp T, 


ow) we, nhy Theoretically, one single equilibrium measure- 
din d In ment in a binary mixture including a determina- 
(3) tion of the temperature is sufficient to calculate 


the two constants, from which the whole equili- 
In case of a two phase-system, eq. (3) is not brium curve can be calculated. 
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An excellent review and critical study of the 
equations mentioned is given by Wont [11]. 

Apart from the limited validity of such a 
calculation owing to the assumption that T and 
p are constant, the equal influence of the vapour- 
liquid compositions and of the temperature on the 
accuracy of the activity coefficient must be seen 
as a drawback. Thus, a possible error in the 
analysis or a departure from equilibrium will be 
enhanced by errors in the temperature. It is 
concluded, therefore, that a testing method 
insensible to errors in the temperature or a 
direct calculation from boiling point-composition 
data is to be preferred. 


EVALUATION OF Vapour-Liquip 
CURVES FROM BOILING 
Data 


Carison and CoLBuRN [1] suggested a procedure 
utilizing the relations 
_ P— Pots and (4) 


P, 2, 

For a given T, 2 diagram y, can be calculated as 
a function of the liquid composition if to a first 
approximation it is assumed that y, = 1, and 
also y, can be calculated assuming y, = 1. By 
plotting these apparent activity coefficients 
against composition the terminal values of y, and 
Ye can be obtained by extrapolation which 
logarithms are the constants in the van Laar 
or Margules equations. 

Orumer, Ricciarp! and TuHakar [4] converted 
the Gibbs-Duhem equation at constant T and p 
to an equation representing 2, y, T at constant 
pressure from which the vapour composition 
could be obtained by stepwise integration. 

A similar relationship was derived by Repuicu 
and Kister [8], it involved however an assump- 
tion whereby the variation of the activity co- 
efficients with temperature at constant composi- 
tion was assumed to be negligible. 

Although the method of Orumer is thermo- 
dynamically exact except for the assumed ideality 
of the vapour-phase the procedure of CarLson 
and CoLBuRN is still to be preferred as this 
method offers a possibility for corrections of 
experimental errors. With Orumer’s method, 


the equilibrium curve becomes, particularly in 
the low and high concentration range, extremely 
sensitive for errors in the 7,2 diagram. 

An exact method will now be derived with 
retention of the possibility for correction of 
incidental errors in the T,2 diagram while 
thermodynamical inconsistencies can be easily 
detected. 

The variation of the activity coefficients with 
temperature at constant composition and pressure 
can be related to the relative partial molal 
enthalpy according to 


diny, 

dT ‘RT? 
where AL,, is the partial molal enthalpy of com- 
ponent 1 in solution minus the enthalpy of the 
pure component at the same temperature. 

As AL = AL, 2, + Alga, where AL is the 
integral heat of mixing per mole of mixture, the 
following relation between the activity coefficients 
of component 1 and 2 can be derived. 


(5) 


a, diny, +a,dlny, = 
ab dT (at pconstant) (6) 
RT? 

By writing eq. (9) as partial derivatives with 
respect to the mole fraction of component 1, a 
modified Duhem equation valid at p constant and 
T variable is obtained. 


diny, din AL MT 
or 
Aln — Aln y, — 
2 1 —2, 1 
1 AL 
er’? (7a) 


The relationship between y, and yg can be 
written as 


in which v = exp —(V, — B,)(P, — p) 
aT » the 
correction term for the non ideality of the vapour 


phase [6, 10}. 
For a given 7,2 diagram y,, can now be 
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calculated for yz = 1 as a function of the liquid 7 
composition and also y, for y, = 1. As an 
example, the binary system methyleyclohexane- | " © 
toluene is given in Fig. 1. The T,2 data of ap 
QuiccLe and Fenske [7] are used for this 7" 
purpose. As the curvature of the calculated curve *7 nD 
is very small, extrapolation to the high and low 
concentration ranges where the data points show y 
026 
02s as 075 10 
3 Mole fraction methyicyclohexane m liquid 
: | \ Fic. 1. Plot of y;, and yg vs. liquid composition for yz = 1 
a20) t and y,; = 1 respectively. 
| 
' | | a large scattering is easily possible. The deviations 
955 —_+ oo in these ranges are caused by inaccuracies in the 
| *;. determination of the 7,2 diagram or departure 
| from equilibrium. 
— . d V1 d z* 
. (8) is 0 
| According to eq. ( rs 7" and o for 
x = 0 and = 1 respectively; thus the inter- 
i r sections of the two curves in Fig. 1 with the 
| 
| | 
| 
4 
: 
acs 
= | 
3 | 
| 
aos 
= 
95 4 ap 020 02% 


Fic. 2. Relation between In y, and In yg with liquid composition as parameter. 
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ordinates at =0 and x= 1 represent the 
maximum endvalues of the activity coefficients. 

In Fig. 2 the logarithms of the smoothed 
pseudo-activity coefficients are used for the 
graphical representation of eq. (8), the parameters 
are taken as ¢ = 0-05 — 0-10... . 0-95. 

It is clear that by stepping over from curve 
m Ar to curve (m+ 1) Ar the accompanying 
change of Aln y,, and Aln y, must satisfy eq. (7a.) 

From the 7, z data and starting with the value 
of In y,, for z = 0 and In y, = 0 obtained from 
the graphical extrapolation in Fig. 1 it is now 
possible, by applying equation (7a) to Fig. 2, to 
draw the entire curve of the activity coefficients 
with composition. 

A simple procedure is to plot the In y,, In y3 
relation on transparent paper, using a graph of 
eq. (7a) with 0-5 Ar — 1-5 Ar... (1 — 0-5 Az), as 
parameter at the same scale, which is put under 
the transparent paper, and then to step off the 
successive curves starting with y, for 2 = 0. 

The graphically extrapolated endvalues of y, 
and y, in Fig. 1 should strictly satisfy the new 
values obtained by eq. (7a) and (8), consequently 
thermodynamic inconsistencies in the boiling 
point - composition diagram can be easily 
detected. 


Aw Exact Test For THE INTERNAL 
Consistency oF Varpour-LiquiIp 
Data 


It is noted by several authors [3, 8] that the 
relation between F*® denoting the excess free 
energy of formation of a mixture over an ideal 
mixture of the same composition, and the relative 
volatility « is eminently suited for testing vapour- 
liquid equilibria. 

However, the authors cited, assumed the 
simple form of the Duhem equation [3] to be 
valid and did not correct for the non ideality of 
the vapour phase. An analogous relation that is 
free of these limitations will now be derived. 

From the definition of the excess free energy it 
follows 


F® = RT (n, In y, + ng In 9) (9) 


By dividing expression (9) by the sum of the 
number of moles and differentiating with respect 


to the obtained mole fraction it takes the form 
dF® diny, d\n 
Introducing the Gibbs-Duhem equation valid at 
p constant 


(10) 


diny, diny,_ AL aT 

az, diy dz 
we see that 

AL dT 


After integration between the limits z = 0 and 
z=1 


f(— fine dx (10a) 


By combination of eq. (10a) and the relation- 
ship for vapour liquid equilibria corrected for the 
nonideality of the vapour phase [6] it follows that: 


An application of this equation to the calecu- 
lation of AF*® from the vapour-liquid equilibria 
data of the system methyleyclohexane-toluene, 
as determined by QuicGLE and Fenske and The 
Central Institute for Physico Chemical Constants, 
Utrecht, Table 1 is illustrated in Fig. 3. 


Table 1. Vapour-liquid equilibria in methyleyclo- 
hexane-toluene at 760 mm Hg. Central Institute for 
Physico Chemical Constants, Utrecht. 


Mole % volatile component in : 
Liquid Vapour Liquid Vapour 
5-00 7-70 55-00 60-80 
10-00 14-60 60-00 65-00 
15-00 21-00 65-00 69-25 
20-00 27-00 70-00 73-55 
25-00 32-60 75-00 77-85 
30-00 42-95 80-00 82-15 
35-00 47-70 85-00 86-45 
40-00 52-20 90-00 90-80 
45-00 56-50 95-00 95-40 
50-00 60-80 


Thermydonamic evaluation of binary vapour-liquid equilibria 


Discussion 
In Fig. 3 In *!, calculated according to the 


method of CarLson and in combination 
with is plotted along with dF*® ‘dx 
against the molal liquid composition. The dotted 
AL . ail from 
RT? dx 
the commonly applied Duhem equation. As 
mentioned before, eq. (3) may hold fairly well in 
a certain part of the diagram and still be in 
serious error at the extremes of the concentration 
At 0-025 mole fraction the value of 


line represents the deviation term 


range. 
mk calculated from Fig. 2 and eq. (7a) is 
dina, 
approximately 0-0125 and 
’ AL aT 
155 
de 
0.02 
| 
r ° 
| | N 
0265 


Mole fraction methylcyclohexane in liquid 
Fic. 3. Plot of In y;/y2 or dF® /dx vs. molal liquid com- 
position for methyleyclohexane-toluene. 


In refers to the method of Cartson and 


only. 
@ method of Cartson and 
{1 method of present author > ats 


A 2, y data of the Central Institute. 
O 2, y data of Quiccre and Fenske. 


Both sides of eq. (3) differ, therefore, by 5-6%, 
when x = 0-025. 

From a comparison of the areas bounded by 
the curve lying above the zero line and the 
curve it follows that by neglecting the 
latter term the overall error in eq. (3) is only 
3%. The deviation of the logarithms of the end 
values calculated from the Margules equation is 
about 

The T, x data of QuiccLE and FeNskr appear 
to be consistent, thermodynamically, and agree 
excellently with the vapour-liquid equilibrium 
data of the Central Institute for Physico Chemical 
Constants. The z, y data of QUiGGLE and Fenske, 
however, are at variance with any equilibrium 
relation. 

It is interesting to note that, on the contrary, 
these authors strongly stressed that conclusions 
as to the validity of equilibrium data should not 
be drawn from temperature, composition data as 
the accuracy of analysis by density or refractive 
index is much higher than those of temperature 
measurements, 
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NOTATION 
A,, Ag = constants in the Margules equation. 
B = virial coefficient. 
FE — excess free energy. 
Jf = fugacity. 
AL = integral heat of mixing per mole of mixture. 


AL, = partial molal enthalpy of component in solution 
minus the enthalpy of the pure component at the 
same temperature. 


n = number of moles. 

P = vapour pressure of pure component. 

p = total pressure. 

R = gas constant. 

T = absolute temperature. 

V = molal volume of component in liquid at tem- 
perature of system. 
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v = correction term for the nonideality of the vapour y = mole fraction of volatile component in vapour. 
phase. y = activity coefficient. 
z = mole fraction of volatile component in liquid. » = chemical potential of a component. 


Subscripts 1 and 2 signify volatile and less volatile com- 
ponent respectively, / and g are abbreviations for liquid - 
and vapour phase. 
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Distillation column for study of individual plate efficiencies 


H. A. C. 
Laboratory of Physics and Meteorology, Agricultural University, Wageningen, The Netherlands 


( Received 15 December 1954) 


Summary —-A laboratory sieve plate column is described for studying the influence of type of 
mixture and operating variables on plate efficiencies. Special care has been given to the possibility 
of regulating and measuring flow of matter and heat to individual trays. 


Résumé— L‘auteur décrit une colonne de laboratoire a plateaux perforés pour étudier l’influence 
du type de mélange et travailler avec des efficacités variables sur les plateaux. II a soigné par- 
ticuli¢rement la possibilité de régler et de mesurer les alimentations et les apports thermiques sur les 


plateaux individuels. 


INTRODUCTION 


A considerable number of methods for predicting 
plate efficiencies in distillation columns has been 
proposed [1, 2, 3, 4, 5]. Many of these are corre- 
lation methods, based upon empirical] relationships, 
others have a theoretical foundation. It is, how- 
ever, necessary to introduce experimental factors 
in the latter cases as well, owing to the complexity 
of the phenomena of mass and heat transfer at 
the conditions prevailing in distillation columns. 

Unless a more fundamental understanding of 
these phenomena becomes available a_ large 
number of experiments continue to be necessary 
if other types of trays are considered, e.g., with 
respect to excessive fouling as may occur with 
some agricultural products, or if mixtures are 
separated which possess properties largely differ- 
ent from those of the mixtures commonly used. 

With the object to obtain data of a fundamental 
character, a laboratory column has been con- 
structed with which the efficiency of individual 
plates can be studied. Special care has been given 
to the possibility of regulating and measuring flow 
of matter and heat for individual plates. Instead 
of using a few trays only as is often done in tray 
efficiency studies, a larger number was deliberately 
installed to enable a comparison of conditions for 
mixtures of different types under nearly the same 
conditions. In the present paper the column is 
described in detail together with an example of 
the results obtained. Complete results have been 
discussed in a previous paper [6]. 
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The column is composed of 5 measuring and 4 
glass sections in alternate position. Each section 
except the uppermost contains a sieve plate. All 
necessary data, viz. liquid flowrate, temperature, 
composition and heat transfer through the wall, 
can be obtained from a measuring section. 

The condenser and reflux divider are mounted 
directly over the top section. The measuring 
sections are insulated by 1 cm of glass blanket, 
the glass sections by means of a double wall. 
Every section is surrounded by a jacket con- 
structed partly of glass, the temperature of which 
can be regulated by means of a set of heating 
coils. Temperature control to within +1°C can 
be achieved. In order to enable a study of 
different types of sieve plates, it was necessary 
to have a free plate surface with a minimum of 
2 rows of holes at a pitch of 15 mm. The limita- 
tions imposed with regard to the minimum down- 
comer and liquid seal dimensions dictated a 
minimum column diameter of 38mm. For the 
present series of runs the hole diameter was held 
constant at 0-75 mm on a 3-O0mm pitch. The 
design of a measuring section is shown in Fig. 1. 
Table 1 summarizes the dimensions of a sieve 
plate. 

For the measurements of the liquid rates on the 
trays a new type of siphon-flowmeter was devel- 
oped which combines small dimensions (7 em) and 
a measuring range of 0-01-lcc/sec with an 
accuracy of about 3%. 
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Table 1. Summary of plate dimensions. 


Inside diameter of column, cm 3-8 
Plate spacing, cm 12-0 
Holes /plate 97 
Hole diameter, cm 0-075 
Drilling pattern Equilaterial triangular 
Distance between hole centres, em 0-3 
Weir height, cm 10 
Inside diameter of downpipe, cm 10 


At closed position of valve V, the liquid flows 
via downpipe D, in the siphon which will fill to 
A and empty. The period of siphoning T com- 
prises the time required to fill the siphon at a 


flowrate of Qcc/sec and a time ¢t’ required for 
emptying; hence T = 6/Q +t’, where b is the 
effective volume of the siphon. A ciné-camera 
served for the measurement of t’. 

The pulsating effect of the siphon on the reflux 
is eliminated to some extent by means of buffer 
tank T. Valve V, can be so adjusted that for a 
siphon content of 4 cc and the dimensions given 
in Fig. 1 the ripple is limited to 8%. These 
fluctuations do not interfere with the operation 
of the column since the siphons are only in the 
circuit during flow measurements. Normally, the 
liquid downflow runs via the downpipe straight to 
the next tray. Even when measuring flowrates 
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Fic. 1. Measuring section of sieve plate column. 
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the disturbances are only local since only one 
siphon is used at a time. 

Pipe E, serves for exhausting the buffer-tank. 
Emptying the siphon has to be accompanied by 
the rapid removal of 4 cc of gas from T,, as other- 
wise an overpressure relative to the siphon inlet 
would develop in the buffer-tank which would 
prevent siphoning. Large values of the equilibrium 
constants, resulting in big temperature differences 
between successive plates, may cause flashing of 
the liquid and, in consequence, a vapour lift in 
the downpipe D,. This is prevented by the exhaust 
pipe E, so the flashed liquid circulates only in the 
circuit formed by E, and T, in which T serves 
as a liquid separator. 

The large diameter of the exhaust pipe is 
required to prevent the formation of liquid 
droplets in it which would hinder rapid exhausting 
and, in consequence, the proper action of the 
siphon. In order to avoid too large holdups, the 
distance between buffer-tank and sieve plate is 
chosen so that at normal throughputs the tank 
drains entirely. It appeared from the temperatures 
recorded that flow-measurements do not influence 
the plate efficiency. 

The results of the flowmeter measurement are, 
naturally, only reproducible by working above the 
so-called ‘* weep-point ” of the sieve trays. In 
order to make possible measurements below the 
weep-point, anti-weeptrays were constructed. 
Their action is based on the principle that the 
pressure in the section below the tray under 
consideration is greater than the pressure at the 
bottom of liquid seal S, so leaking liquid will be 
driven back through the drainage holes via the 
liquid seal to the tray. Since the liquid runs as 
a film along the surface of the tray to the edge 
and there drips off, a narrow annular model as 
shown in the figure was sufficient. 

Sample taps specially designed for withdrawal 
of small samples (0-2 to 0-4 cc) are provided on all 
measuring sections and on two of the glass ones. 
The vapour sample inlet of section 6 is located 
about 05cm below the anti-weep tray and 
exactly in the centre of it. The liquid remaining 
in the sample line is blown back into the column 
just before sampling. 

Every section is provided with a thermo- 


element. The temperatures of the vapour as well 
as the temperature differences between the section 
and its heating jacket were measured by means 
of a specially designed circuit. In order to avoid 
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Fic. 2. Flow sheet of laboratory distillation unit. B S.F.— 
buffered siphon flowmeter; F - flowmeter; S — sample 
tap ; T thermocouple ; R — rotameter. 
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an intermittent distillation caused by boiling 
retardation and non-mixing in the reboiler, a 
mixer was installed in it. 


FLow SHEET 


A flow sheet for the operations is given in Fig. 2. 
Overhead vapour is wholly condensed. For runs 
under finite reflux a part flows by gravity via 
pilot valve and flowmeter to the reboiler. For 
runs with total reflux for one component only the 
distillate is removed as product. Total reflux for 
one component only is achieved by supplying one 
pure component from a constant head tank via 
the measuring and regulating system to the 
reboiler, the other component from a similar 
system to the middle of the condenser, in order 
to avoid upsetting of the heat balance and the 
plate efficiency, and removing as much distillate 
as amounts to the total supply of pure com- 
ponents. The upper section does not contain a 
sieve plate and the reflux from the condenser 
flows directly through downpipe or siphonmeter. 
The supply of pure components and removal of 
distillate is controlled by means of completely 
buffered siphon-flowmeters. The flowrate of 
cooling water in the condenser is controlled with 


the aid of a rotameter. A continuous record of 
the vapour temperatures of all sections is obtained 
by means of a 12-point recorder. 

A run at total reflux of the volatile component 
only is illustrated in Fig. 3. 


Z 


02 Q4 06 ae 10 
=~ Mole fraction methyicyclohexane in liquid 
Fic. 3. MeCabe-Thiele diagram of methylceyclohexane- 


toluene at partial reflux of the volatile comp. only, 
L/V = 0-84. 
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Rectification of multicomponent mixtures 


A. M. Esnaya 
The Coca-Cola Export Corporation, 515 Madison Avenue, New York, N.Y. 


(Received November 1954) 


Summary—On the assumption that the relative volatilities, «, and the molal overflow are 
constant in a given column section, a comparatively simple and rapid method of performing 
multicomponent fractionation calculations is developed. The proposed method enables one 
to calculate by means of tabulated values of the functions E, and Ey directly the composition 
on the nth plate without going through calculations of the compositions on the intermediate 
plates. This is done by means of the simple equations : 
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”) (in stripping section) 


values may be obtained from suitable tables or plots. 

It is believed that the proposed method offers considerable simplicification in multicomponent 
rectification calculations. “he final equations are tested and their use demonstrated in two 
examples. 


Résumé—En supposant que les volatilités relatives, «, et le flux moléculaire sont constants 
dans une section donnée de la colonne, l'auteur expose une méthode comparativement simple 
et rapide pour effectuer les calculs dans la séparation fractionnée d'un mélange a plusieurs 
composants. La méthode proposée pern.et, & partir de tables des valeurs des fonctions E, et Eg, 
de calculer directement la composition sur le niéme plateau sans passer par le calcul des com- 
positions sur les plateaux intermédiaires. Ce qui se fait par les équations simples : 


+ p, (in rectitying section) 


= Ey + (in stripping section) 

E, et E, sont des fonctions du taux d’équilibre du reflux et du nombre de plateaux. Leurs 
valeurs peuvent étre obtenues a partir de tables ou de courbes appropriées. 

L’auteur pense que la méthode proposée présente une simplification considérable dans les 
calculs de la rectification de mélanges & plusieurs composants. II essaie les équations finales 
et démontre leur utilisation dans deux exemples. 


The classical method of solving multicomponent 
rectification problems is that of Lewis and 
MaTueson [8]. This method consists of plate-by- 
plate calculations starting with the still pot and 
continuing up the column until a composition 
approaching that of the feed is obtained. In a 
similar way, starting with the composition of the 
distillate, one proceeds down the column until a 


composition close to that of the feed is reached. 
Readjustment of the reflux ratio or of one of the 
starting compositions (waste or distillate) may be 
required to make the compositions “ mesh ” on 
the feed plate. This method (as well as its graphical 
equivalent, the method of Lewis and Cope [2]) 
has the advantage that it may be carried to any 
degree of accuracy and may be easily combined 
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with enthalpy balances in the cases in which the 
assumption of equal molal overflow is not justi- 
fied. The only disadvantage of these classical 
methods is the tediousness of the calculations 
involved: a trial and error procedure in the 
calculation of the composition of the vapour 
rising from each plate. 

Because of the practical importance of the 
problem numerous workers in the field have 
proposed methods which are intended cither to 
simplify the necessary calculations, or to give an 
answer to some specific aspect of the general 
problem. Important recent contributions have 
been that of Murpocu and co-workers [4, 5, 6] 
who have presented a mathematica! analysis of 
the general problem, and that of Unperwoop 
(8, 9, 10, 11, 12] who has developed a rigorous 
solution for the composition on the n' plate 
without going through plate-by-plate calculations 
for all intermediary plates. The method involves 
the solution of a system of equations, each system 
having as many equations and rootes as the 
number of components involved. The amount of 
numerical calculations required is thus con- 
siderable (see ref. 8). 

The present method allows the calculation of 
the composition of the n™ plate directly in a 
convenient way by means of tabulated functions. 

One must consider first the following simplified 
problem : 

A mixture of z components enters the feed plate 
and the composition of the mixture is known to 
be ty, ty, (Note: the first subscript 
refers to the component and runs from a to 2; 
the second subscript refers to the number of the 
plate). The composition of the distillate is r,,. 


Tyg and the composition of the waste 


R> Ri ic. there are no “pinched in” 


regions in the section. The problem is to evaluate 
analytically the number of plates in the rectifying 
and stripping sections of the column. 

Let us consider any one single component of 
the mixture. We can drop the subscripts a, b, 
. . . z and write for the equation of the equili- 
brium curve : 


y, = (1) 
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and for the operating line : 


= — 4 (2) 


where p is the slope of the operating line in the 
section of the column under consideration, and 
— a is the intercept of the operating line with 
the 2 = 0 axis. 

Using equations (1) and (2), one can show that 
the composition on the nth plate is given by : 


p” 


Equation (3) is basically the same as the one 
derived by Epsisrer [1]. Implicit in it is the 
usual assumption of constant molal overflow 
since p has been considered as a constant in the 
section under consideration. The equilibrium 
ratios K, as shown in equation (3) are functions 
of temperature {at a given pressure) and of the 
composition. In order to obtain a convenient 
analytical solution, it is necessary to assume 
ideality of the liquid phase. This assumption 
makes the equilibrium ratios K, independent of 
composition. 

In order to evaluate anenytenny Ge products 


of equilibrium ratios of the type 7 '*K, entering 


into (8), it is necesssary to etablish some sort of 
relationship among the various K’s. Before trying 
to determine the simplest suitable form of this 
relationship it should be noted that what deter- 
mines the composition of the vapour in equili- 
brium with the liquid phase on an ideal plate are 
not the absolute values of the K’s of the various 
components, but rather their ratios to the K of 
one component selected at random. These ratios, 
known as relative volatilities, 2, are the sole 
factors in determining the composition of the 
vapour. Thus, one can write all the pertinent 
rectification equations in terms of « alone without 
the K values entering into the picture. In accord- 
ance with the practice accepted by recent con- 
tributors in the field, the relative volatilities, «, 
will be assumed constant in the section of the 
column under consideration. 


From the above it becomes clear that the 
relationship connecting the K’s of the same 
component on successive plates must meet only 
the following three requirements : 

(a) It must give on every plate such values of 
K, that the ratios of the K’s for the various 
components are the correct ones. 

(6) In view of the assumption of constant 
relative volatility the ratios of the K’s must 
remain constant on successive plates for any 
given pair of components. 

(c) It must contain an adjustable parameter 
which will make it possible to satisfy the criterion 

= 


2 x; = 1-000, 
a 


A relationship which meets the above require- 
ments is : 


K,=K,b‘ (i =1,2,8,...n—1) (4) 
provided b is the same for all components because: 


K.. = Keo by = Keo = 4, constant (5) 


which means that 6, = 6, = 6 for all i + 0. 

Equation (4) should therefore be looked upon 
only as a convenient mathematical device for 
meeting the requirement of correct constant value 
of « on every plate. How closely this equation 
approximates the true curve of K vs. i in the 
column is immaterial, and the values of K, given 
by it are not and do not have to be the true 
ones on the successive plates. The proposed 
method, therefore, contains no other assumption 
than that of constant molal overflow and constant 
relative volatility and gives exact results when- 
ever the molal overflow and the relative volati- 
lities are constant. 

With the relationship expressed in (4) the value 


of the products a K;, is equal to : 
4 


n-1 -Zi 
q 


K,= (6) 


and equation (3) may be rewritten as : 
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Kop 


p” 


n—2 —mn—s) + 


ar 


Referring all relative volatilities to component 
a, we have : 


a 
my +S {1+ 


(7) 
+ 


K, = K,; K, = K,%,; 
K, = = K,4,,; 


and (6) when written for component z becomes : 


nm K,, = (8) 
i=q 


One can now obtain the final form of equation 
(7) which is : 


nin—1) 


n—2) (n— 9) 
( 
4 
where : 
Ka 
Q= 
Pp 


If one denotes by E, the expression Q" b~~ 
and by E, the expression in the large brackets, 
one can rewrite (9) as follows : 


The values of a and p are : 
in stripping section: a, = 2, (p, — 1) 


(slope of operating line 
in stripping section) 
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in rectifying section: a, = 2, (p, — 1) 
V 
(slope of operating line 
in rectifying section) 


Thus, if one denotes by m the number of plates 
in the stripping section (counted from the re- 
boiler) and by n the number of plates in the 
rectifying section (counted from the feed plate,) 
the following equations apply : 


2, = 2, EB, 4 (P, — 1) E, (11) 
s 
x, =2,E, + (12) 
If equations (11) and (12) are applied to a sub- 
section not starting at the bottom of the stripping 
or rectifying sections, their form becomes : 


= % E, + E, (13) 
=a E, E, (14) 


where 2, denotes the first plate of the 
subsection and m and n give the num- 
ber of plates in the subsection, i.c. 
above 2,. 


One should consider now in greater detail the 
functions EZ, and E,. Both of these functions 
depend on Q, b, and n (or m). Table 1 gives the 
forms of E, and E, for m (or n) from 1 to 10. 

Both E, and E, can be calculated and plotted 


as functions of Q with b as a parameter. Such a 
plot for n = 5 is shown in Fig. 1 on a greatly 
reduced scale. 
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Tabulated values of E, for Q = 0-5, 1-0, 1-5, 
2-0, 2-5, 3-0, 3-5 and 4-0; b = 1-0, 1-01, 1-02, 1-08, 
1-04, 1-05, 1-06, 1-07, 1-08, 1-09, 1-10, 1:15 and 
190; « = 2,8,.... 10 are given in Appendix I. 

The method of determining the composition on 
plate m (or ») is relatively simple. Noting that 
K,, is the equilibrium ratio of the component to 
which all the relative volatilities are referred and 


Table 1. E, and E, for n or m equal to 1, 2, .. . 10. 


norm E, 
1 Q 1 
2 er 
3 1+@? +. 
4 14+@*%* +. 
5 1+Qo* +@o7 + + Qn? 
7 1+@Q@* 
1+ Qo” 
9 1+@* 
10 1+@* 
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at the temperature of the bottom plate of the 
section or sub-section considered, all the Q’s for 
the various components are computed. Thus, if a 
is the component to which all the relative volati- 
lities are referred, one has : 


Q, p Q, ™ 


Depending on whether the stripping or rectify- 
ing section or a subsection in either is being con- 
sidered one of equations (11) to (14) inclusive is 
written for every component. Thus if the stripping 
section is being considered the equations written 
are : 


= E, (Qa, b, m) + 


= E, (Qe, b, m) +. (Q., b, m) 


The values of p and 6 are the same for all com- 
ponents. 

Once equations (14) are written, one selects a 
value of 6 at random and determines from the 
plots of E, and E, the values of these two functions 
for each component and for the particular value 
of m (orn). These values of E, and E, are used in 
equations (15) and Dims Calculated. 
If these satisfy the condition 


Ze, = 1-090 (16) 


then the value of b used was the right one and the 
concentrations 2,,,, Zym_, ete., are the correct 
ones. If (16) is not satisfied, a new value of 6 is 
selected and E, and E, corresponding to this new 
value are read off the charts. The calculation 
may thus have to be repeated several times before 
the value of b which satisfies (16) is obtained. The 
process is a very rapid one, however, since on 
every trial the only variable terms are E, and E, 
while the coefficients in front of them remain 
the same. 

The present paper gives values of E, for n from 
1 to 10. This limits at present the number of 
plates which can be considered at a time to 10. It 
is possible, however, to extend in the future the 
range of known values of E, and thus make 


possible the calculation of 20, 30 or more plates 
in a single step. In general the values of p and « 
cannot be considered constant in sections con- 
taining more than ten plates except, perhaps, in 
some very difficult separations. It is advisable, 
therefore, for greater accuracy to calculate the 
column in subsections of ten plates, check the 
temperature and reflux conditions prevailing on 
that plate and use the new values of p and « thus 
determined in the calculation of the next sub- 
section. It should be noted that in order to 
calculate K,, at the plate on the bottom of the 
following subsection (equal to K,,, at the top 
plate of the preceding subsection) one should 
not use the formula given in (4), but rather 
calculate the temperature on the plate and use 
K,,, corresponding to this temperature. The 
reason for this is that, as previously explained, 
equation (4) does not necessarily give the correct 
values of the K;’s but only the correct values 
of a. 

The value of b is not necessarily the same for 
two consecutive subsections. The component to 
which all relative volatilities are referred should 
be one of intermediate volatility so that it is 
present in significant amounts in both ends of the 
section (or subsection considered). 

One of the basic difficulties in multicomponent 
rectification calculations is the impossibility to 
fix all terminal composition once the feed com- 
position and the reflux ratio have been fixed. (For 
detailed discussion of this aspect of multicom- 
ponent rectification see ref. 7). This difficulty of 
uncertain terminal composition arises from phase 
rule considerations and no method can ever hope 
to overcome them. All present methods handle 
them by making reasonable assumptions as to the 
concentrations in the terminal streams, and 
checking these assumptions against the results 
obtained. The method proposed in this paper 
should be used in the same manner. 


Minimum number of plates 
For this special case p = 1 and equations (11) 
and (12) become : 

= 2, E, = 7, b- (17) 
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In determining the minimum number of plates 
required in order to reach a certain composition 
both z,, and 2, (or 2, or a,) are given, and (17) 
may be rewritten as : 


b= in stripping section 
in rectifying section 
Equations of the type of (18) may be written for 


every component and since 6 is the same for all 
components, one obtains : 


I II 


It can be shown that in the case of only two 
components equation (19) is identical to the 
Fenske equation. 

One can take I and II and solve directly and 
easily for the unknown nx. Then in a similar way 
take II and III, IIT and IV, ete. If the values of 
n thus obtained differ considerably, it shows that 
it is impossible to obtain at total reflux the con- 
centrations 2. Starting with the 
concentrations 14. 7 A readjustment 
of one of the terminal streams is then made and 
the calculation repeated at least until » for the 
important pairs of components is the same. Since 
the number of plates at minimum reflux is directly 
determinable from the desired separation of the 
keys, the important pairs of components are 
clearly the key components. The illustrative 
example shows the manner in which the minimum 
number of plates is obtained, or conversely the 
manner in which the composition on the n™ plate 
at total reflux is determined. 


Minimum reflux ratio 


Determination of the minimum reflux ratio in 
multi-component rectification is one of the most 
difficult aspects of the general problem. Equations 
(11) and (12) can presumably be used with 
successively decreasing values of p, or increasing 
values of p, until a value of p, (or p,) is reached 
which satisfies the criterion for minimum reflux. 
Since p, however, enters in both E, and E, as 
well as in the coefficient in front of E, (Eq. 11 
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and 12), the procedure is slow, and the present 
methed is not recommended specially in the 
presence of other excellent methods for obtaining 
the minimum reflux ratio [4, 5). 


Example No.1 


As a first example we shall consider the tar 
acid fractionation problem solved by plate-to- 
plate calculations by Rorinson and 
[7]. p. 236, 

The composition of the waste, feed, and dis- 
tillate streams are as follows : 


Waste Feed Distillate 
0-0524 0-350 09530 
0-2020 0-150 0-0455 
0-4470 0-300 0-0015 
0-22.40 0-150 
r 0-0750 0-050 


The slope of the operating line in the stripping 
section is 1-184 and in the rectifying section 
0-9091, The temperature in the reboiler is 160°C 
and at this temperature, the relative volatilities 
are 


125; 4, 100; a, = 0-70; a, 0-44; 
x4 — 0-087. K, at 160°C is 1-46. 


Since the purpose of this example is to compare 
the results of the present method with the 
rigorous stepwise method, the procedure used by 
Rostxson and will be followed closely. 
Thus, since these two authors change the values 
of x on the 8th plate above the reboiler, a sub- 
section starting with the still pet and including 
the 8th plate above it will be calculated first. 


= = 1231 125 = 154 
& = 1231 x 100 = 1-231 
® = 1231 X 0-70 = 0-862 
Qa = 1231 K 0-44 = 0-542 
= 1-231 0-087 = 0-107 


Equations (11) for all component on the 8th 
plate becomes : 
0-0524 E, (1-54, b, 
= 0-202 b 
0-447 


+ 0008153 Ey (1-54, 8) 
0-03143 Ey (1-231, b, 8) 
006047 E, (0-862, b, 8) 
0-03485 Ey (0-542, b, 8) 
0-01167 Ey (0-107, b, 8) 
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e 
By trial and error it is found that J 2, = 1-00 


a 
when b = 1-066. Substituting the values of E, 
and E, at b = 1-066 in the above equations, one 
obtains the composition in the 8th plate : 


Zag = 0-0524 X 5-250 + 0-008153 R. and G. 
X 18-82 = 0-2751 + 0-1084 = 0-3835 0-375 

pg = 0-202 X 0-8806 + 0-03143 
X 5:59 = 0-1779 + 0-1757 = 0-3536 0-355 

= 0-447 X 0-0508 + 0-06947 
x 246 = 0-0227 + 0-1711 = 0-1938 0-204 
= 0-224 X O + 0-08485 x 1:59 0-0552 0-056 
= 0-075 X 0 + 001167 X 1-16 = 0-0135 0-013 


a2 = 09996 x = 1-003 


The temperature on the 8th plate may be found 
to be 120°C and Ropinson and GILLILAND use 
the following values of the relative volatilities for 
the next subsection : 


Xp, = 1:26; a, = 1:00; a, = 0-675; ay = 0-392; 
%& = 0-087. K,, (at 120°C) = 1-013. 


With these values of « and K,, the Q’s for each 
component are calculated and one can write for 
the composition of the 5th plate of the sub- 
section which starts at the 8th plate above the 
reboiler (i.e. the 13th plate counting from the 
reboiler) the following set of equations : 

= 03835 E, (1-078, b, 5) + 0-008153 Ey (1-078, b, 5) 
Zps5 = 0°3536 E, (0-8553, b, 5) + 0-03143 Ey (08553, b, 5) 
= 01938 EB, (0-5773, b, 5) + 0-06947 (0-5773, b, 5) 
= 00552 E, (0-3353, b, 5) + 0-03485 Ey (03353, b, 5) 
2,5 = 0-0135 E, (0-0744, b, 5) + 0-01167 Ey (0-0744, b, 5) 


After several trials it is found that the value of 
e 
b which satisfies Yr, = 1-0 is b = 1-015. With 


a 
this value of b, the functions E, and E, are read 
off the charts, substituted in the above equations 
and the composition on the 5th plate of the 
subsection obtained : 


Robinson and Gilliland 


0-5242 0-525 
= 0-2527 0-252 
243 = 01601 0-160 
Tyg = 0-0528 0-0514 
= 0-0126 0-0126 
-= 1-0024 Zr = 1-0010 


Because of the fact that the ratio of the con- 
centrations of componnet a to component 6 in 
the liquid on the 18th plate is essentially that of 
the feed, Rosinson and consider the 
13th plate as the feed plate. The same will be 
done in the present example. In the rectifying 
section p, = 0-9091 and K, ,, = 0-9605. 


= =1-0565X 1-26 =1-8312 
= 1-0565=1-00 =1-0565 
Q. 1.0565 x 0-675 =0-7181 
= 1-0565 x 0-392 =0-4141 
= 1-0565 x 0-087 —0-0019 


Zan = 0°5242 E, (1-8312, b, n) — 0-0958 (1-812, b, n) 
Zp, = 0°2527 E, (1-0565, b, n) — 0-00455 E, (1-0565, b, n) 
Zen = 01601 E, (0-7131, b, n) — 0-00015 Ey (0-711, b, n) 
Zam = 00528 E, (0-4141,6,n) —0 Ezy (0-4141, b, n) 
Zn, = 00126 E, (0-0919, b,n) —0 X Ezy (0-0919, b, n) 


To calculate the composition on the 10th plate 
above the feed plate (23rd above the reboiler) 
one uses the values of the functions E, and E, 


e 
at n = 10. The value of b which satisfies 2 zy = 


a 
1-000 is b = 1-0257. The equations for this plate 
then become : 


* 10 = 0°5242 X 5-578 R. and G. 
— 0-0953 X 21-38 = 0-8865 0-896 

Zp, 19 = 02527 X 0-5532 

— 0-00455 X 6-871 = 0-1084 0-102 
= 01601 X 0-01085 

— 0-00015 X 2-483 = 0-0013 0-0007 
= 00528 x 0-00004 = 0-0000 0-000 
2, 9 = 0-0126 x 0 = 0-0000 0-0000 


Zz = 09962 = 0-9987 


Next one can calculate the subsection included 
between the 10th and 13th plate above the feed 
(23rd to 26th plate above the reboiler.) The new 
values of Q are : 


Q= =0-8922 x 1-00 ~0-8922 
Q.= =0-8922 x 0-675 =0-6022 


The composition on the 13th plate above the 
feed is given by . 
7, = 0-8865 E, (1-1242, b, 3) — 0-0953 E, (1-1242, b, 3) 


2, = 01084 E, (0-8922, b, 3) — 0-00455 E, (0-8922, b, 3) 
2, = 0-0013 E, (0-6022, b, 3) — 0-00015 E, (0-6022, b, 3) 


ALLEN M. Esnaya 


The value of 6 which gives J'z, — 1-000 is 
= 3-49 and with this ene obtains : 


2q = 0-8865 X 1-4208 — 0-0953 R. and G. 
X 83-8880 = 1-2595 — 0-3229 = 0-9366 0-942 

2» = 0-1084 x 0-7102 — 0-00455 
X 26882 = 0-0770 — 0-0122 = 0-0648 0-058 

z, = 0-0013 x 0-2184 — 0-00015 
X 19648 = 0-00029 — 0-00029 = 0-0000 0-000 


= 10014 Ir = 1000 


The vapour rising from the 26th plate upon 
condensation gives the desired distillate com- 
position, and the column therefore has 25 plates, 
18 in each rectifying and stripping section. The 
answer in the above example was obtained in 
four trial-and-error steps instead of the 26 such 
steps of the rigorous step-wise calculation. The 
considerable saving in time is thus obvious. The 
relationship between the actual values of the 


+0100 \ 
« 


volatility ratio and those given by K,; = Ky x 
1-066" is shown on Fig. 2. 


Example No. 2. Minimum number of Plates 


Consider the following five components under 
total reflux : 


Component x, (assumed) 
a 0-0524 1-25 0-585 
b 0-202 1-00 0-376 
c 0-447 0-70 0-021 
d 0-224 0-44 0-018 
e 0-075 0-087 0-000 


What is the minimum number of plates required 
to obtain the composition 2, (assumed) ? 
From equation (19) we have : 


0-585 3-376 
I 
2. x 0-70") ( 
0-021 0-018 
IV 

From I and II: n=8 
From II and III : n= 10-4 
From IIland IV: n< 0 


Clearly the assumed concentrations cannot be 
obtained under conditions of total reflux and the 
composition on the n™ plate should be readjusted. 

If we consider components a and b as the 
important ones and do not wish to change 2,, 
and 2,,, we can calculate readily the change 
necessary in the assumed values of x, and x, in 
order to obtain a composition compatible with 
the starting composition under total reflux. Thus 
solving IT and III together with n — 8 we obtain : 


(Sere x 1-008) 0-70") 
0.376 
From II and IV: 
0-876 Tan 
Tan 0-0006 


The results are summarized and compared with 
these obtained by stepwise calculations follows : 


VOL. 
4 
1955 
ACTUAL VALUES OF 
K = Ky x 
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8 
Component (stepwise calculation) (present method) 
a 0-585 0-585 
b 0-376 0-376 
c 0-048 0-048 
d 0-000 0-0006 
e 0-000 0-000 


NOMENCLATURE 
a4 = intercept of operating line with « = 0 axis 
4, = intercept of operating line in stripping section with 
x = 0 axis, a, = 2, (p, — 1) 
a, = intercept of operating line in rectifying section with 
a = Oline,a, = —2g(1 — p,) 
b = a constant in the equation K, = Ky 6"; « para- 
meter in the functions E, and EB, 


i = running index for number of plate — 0.1.2... 

m = number of plate in stripping section starting with 
reboiler 

n = number of plate in rectifying section starting with 
feed plate 


p = slope of operating line 
Pp, = slope of operating line in stripping section ; 


L 


P, = slope of operating line in rectifying section ; 
L L 
# = concentration of a component in the liquid phase, 
mole fraction 
z, = concentration of a component in the waste, mole 
fraction 
ay = concentration of « component in the feed, mole 
fraction 
7g = concentration of a component in the distillate, mole 
fraction 
x, = concentration of a component in the liquid on the 
plate 
% = concentration of a component in the liquid on bot- 
tom plate of a subsection 
y = concentration of a component in the vapour phase 
D = moles of distillate 
E, and E, = functions defined in equation (9) 
K = equilibrium ratio = y/x 
L = moles of liquid overflow (1 moles of liquid overflow 
in stripping section) 
P 
}’ = moles of vapour going through column 
W = moles of waste 
2%) = relative volatility of component a referred to com- 


ponent b ; ab = 
b 
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Appendix I 


TasuLaTeD VaLues or E, (Q, 6, n) 


for: n = 2, 3, 4, 5, 6,7, 8, 9, 10 


Q = 0-50, 1-00, 1-50, 2-00, 2-50, 3-00, 3-50, 4-00 
b = 1-00, 1-01, 1-02, 1-03, 1-04, 1-05, 1-06, 1-07, 1-08, 1-09, 1-10, 1-15, 1-20 


Note: E, (0,6, n) = 1-000 
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n=2 
b | Q= 050 | 100 Q= 150 | 200  Q= 250 | 300 | Q= 350 | Q 400 
1-00 15000 20000 2500 8-000 4-000 4-500 5-000 
101 | 1495 | 1990 | 2485 | 2080 8475 3-970 4-465 4.960 
102 1-490 1980000 2471 | 2061 3-451 3-941 4431 4-922 
1-08 | 1-485 1-971 2456 2062 | 3-427 3-913 4-308 4884 
104 148162 2-923 3404 4.366 4846 
1.05 | 2620 2-005 3-381 3-857 4333 4510 
1-06 2-887 3-358 3-830 4-302 4774 
107 1935 | 2402 2-869 3-336 3-804 4271 +738 
1-08 1463 | 1-926 2-388 2-852 3-315 3-778 4241 4704 
1-09 1459 | 2-376 2-835 3-204 3-752 4211 4-670 
1-10 | 1-455 | 1908 264 2818 3-273 8-727 4182 4636 
1-15 1-435 1870-2807 2-739 174 3-609 +044 4-478 
1-20 2-667 3-083 3-500 3-917 4-333 
ne8 
b | 100 Q= 1H 200 | = 250 @= 300 | 3:50 Q = 400 
1-00 | 1-750 3-000 47 000 | 9-750 13-00 16-75 21-00 
101 | 2-951 4654 6-848 9-517 12-68 16-32 20-45 
102 1-716 2.904 4-562 6 692 9-202 12-36 15-91 19-92 
1038 1-700 2.858 4473 6546 0-076 12-06 15-50 19-41 
104 2814 4-387 6-405 8-878 15-13 18-92 
1.05 2-771 +304 6-269 | 8 666 11-50 14-76 18-45 
1.06 1-655 2-730 224 8666-138 8472 11-23 14-40 17-99 
107 1-641 2-690 4147 8-286 | 10-97 14-06 17-55 
1-08 627 2-651 4072 58900 8105 10-72 13-72 17-13 
1-09 1616 2614 4-000 5772 | 7090 «(10-47 13-40 16-72 
110 1-601 2-578 3-930 5-658 7762 «10-24 13-10 16-32 
1s 3-614 | 7000 11-70 1454 
1-20 2-273 4708 | 6-353 | 8-292 | 10-52 13-04 
| 
ne4 
| | 
b 10 BH 200 | so | Q = 400 
1.00 (1875 4-000 8-125 1500 25-87 40-00 | 59-63 85-00 
101 3864 7-776 1428 | 24-09 S791 80-40 
102 3-736 7 19610 2289 76-09 
1-03 1-778 3-615 6140 1298 | | 72-06 
104 | 17490 6-851 1230 (82400 4807 68-29 
1-05 1-721 3-394 6377 | #1283 | 1971 | 3079 6475 
1-06 1695 3-292 6320 1131 1878 | 614s 
107 1-670 3-196 6-078 10-82 17-91 27-86 
1-08 1646 8-105 5849 (10:85 17-08 26-52 | 
1.09 1623 8.018 5633 9914 1630 2526 | 8723 | 
1-20 143200 2816 8003 1019 | 31-18 


oF 
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neS 
| 
b Q = 0-50 Qu 100 | Q = 200 @-250 | @-s00 | Q = 400 
1-00 1-937 5-000 | 13-19 3100 64-48 2098 341-0 
1-01 1885 4713 | 12-21 28-45 | 58-88 nos | 1909 310-0 
1-02 1835 4452 11-32 2615 3387 102-0 173-8 282-2 
1-03 1790 4212 1052 24-06 49-33 93-75 | 158-6 257-1 
104 | 4-747 3-993 9-783 22:18 | 45-25 86-28 | 144-8 284-5 
1-05 | 3-792 9-116 20-47 41-55 76-99 1824 2141 
1-06 1671) 38007 8-508 18-91 38-19 73-28 195-8 
107 1637 § 3438 7-955 17-50 35-16 6764 | 110-9 179-0 
1-08 1-605 3-282 7-449 16-22 32-39 62-49 101-7 163-9 
1-09 1-575 3-138 6-985 1505 | 20-88 52-77 93-32 150-2 
1-10 1-547 3006 6-562 13-98 | 27-60 48-51 85-74 137-8 
1-15 1-431 2-479 4914 9876 | 18-87 33-79 56-93 90-77 
1-20 1-345 2-117 3-823 7-215 | 13-28 23-27 38-66 61-14 
| 
n= 6 
b Q= 050 | Q@=100 | 150 @=200 | @=250 | Q = 400 
1-00 1-969 6-000 20-78 63-00 162-1 3640 8 § 7349 1,365 
1-01 1-897 5485 | 18-42 55-14 | 141-1 3158 | 6366 | 1,181 
1-02 1-831 5032 1638 | 128-0 2745 | 5522 | 1,023 
1-03 1-772 4633 | 1460 #=‘| 4252 | 107-4 2389 479-8 888-0 
1-04 1-718 4-382 1306 | 37-46 92-99 208-4 417-6 772-0 
1-05 1-669 3-970 | 11-710 || 88-07 82-38 182-0 | 363-9 671-9 
1-06 1-624 3695 9536 29-26 72-34 159-2 317-6 585-7 
107 1-584 | 3-452 9-509 25:96 6367 18950 SS 
1-08 1-547 3-234 8605 23-08 | 56-12 1225 | 234 | 4872 
1-09 1-512 3-040 7-810 20-56 49-55 107-6 2132 | 
1-10 «(18-36 43-85 9479 | 1873 | 38422 
1-15 1-356 2232 4065 | 10:77 | 23-46 5140 99-99 181-5 
1-20 | 1-270 1-851 3-305 6-800 | 14:35 | 29-06 | 55-87 | 99-30 
b Q= 050 | Q=100 | | Q=200 Q= 250 Q= 300 | Q=350 | Q= 400 
1-00 1-984 7-00) 32-17 127-0 | 406-2 1,093 2,573 5,461 
1-01 1-893 6-167 27-04 1049 = 883-2 893.6 2,100 4,451 
1-02 1-813 5-468 22-82 86-90 274-0 782-2 1,717 3,636 
1-08 1-742 4-880 19-35 72-30 225-8 601-3 1,407 2,976 
1-04 1-679 4-384 16-49 60-22 186-7 495-1 1,156 2,442 
1-05 1-617 3-981 4 50-36 154-7 408-4 951-4 2,006 
1-06 1572 | 3-605 12-14 42-26 131-5 337-6 784-6 1,651 
1-07 1-528 | 3-300 10-50 35-60 107-0 279-9 648-7 1,364 
1-08 1-487 3-038 9-134 30-09 89-42 232-5 587-5 1,127 
1-09 } 9-451 2-812 7-985 25-52 74-86 193-5 446-0 934-6 
1-10 1-418 2618 7-022 21-73 62-87 161-5 370-9 775-9 
1-15 1203 1-965 | 4-024 10-35 27-43 67-65 152-2 3148 
1-20 1-213 | 1-620 2-660 5-553 13-01 30-18 65-84 133-9 


OL. 
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n=8 
| | 
b | 080 | @=100 Q=200 | 250 | 300 Q= 350 Q= 400 
1-00 255-00 10166 3,280-0 9,006-9 21,845 
101 1-883 196-69 778-05 | 2.5015 6856-1 | 16,609 
1-02 1789 760 30-801 152-31 597-36 5,233-3 -12, 662 
1-03 1708 4068) 24-506 118-42 460-09 | 1471-6 4.005-4 9,678-7 
1-04 1638 4-332 | 19-794 92.520 355-74 | 1120-6 3,075-9 7422-4 
1-05 16-048 72-579 | 275-86 871-70 | 2,366-7 5,705-3 
1-06 57-2080 21464 1827-6 4,396-9 
107 1476 10-813 45-338 | 16770 52396 | 14152 8,399-2 
1-08 36-104 | 131-42 407-94 2,633-6 
1-09 1397 7-552 28-912 103-37 854-70 2,045-4 
1-10 | 23-329 81-726 249-84 667-76 | 1504-7 
1-13 3-271 87922 26-813 277-408 201-62 475-12 
1-20 1169 1-852 2014 41015 | 10-083 | 26-286 65-373 150-59 
| | | 
n= 
| | | 
| 
1-00 748900 SO 2,542 31,530 87,370 
101 1869 7-286 5479 8643 | 1,797 «6,931 22,160 61,350 
1-02 | | 5-016 043 210 1,275 15,630 43,220 
1-03 1674 #923 30-12 1880) 909-1 BATT 11,070 30,560 
1-04 1599 4171 22-84 13977 6386 2,515 7,997 22,070 
1-05 1-534 3-588 17-29 99-24 4678 «(1,971 5,609 14,740 
1.06 1-478 3-132 13-34 7278 4,014 11,030 
1-07 1-429 2-778 10-44 53-79 | 2451 9163 2,885 7,918 
1-08 1-387 2-498 8-290 wuz 1785) 2,078 5,692 
1.09 1-350 2-274 6685 3002 | 130-7 1,502 4,106 
1-10 1318 2-093 5483 | 96-27 | 350-5 1,001 2,975 
115 120821368 | 76-84 231-4 621-5 
1-20 1-136 | 1338 1-738 | 2-910 6-873 19-37 54-29 141-2 
| | | | 
n= 
b @=030 | @=100 | | | 2% Q=300 Q=350 | Q=400 
1-00 | 1998 1,023 6,358 29,520 110,300 | 349,400 
101 667-2 4,109 19,010 70,920 | 224,400 
1-02 | 5-951 51-74 437-8 2.669 12,300 45,790 144,700 
1-03 | 35-63 289-1 1,743 7,904 29,690 93,680 
1-04 1-562 3-927 24,92 192-4 1,145 5,215 19,360 60,990 
1.05 1-494 3-310 17-772 1200 T3547 38,426 12,660 89,830 
1.06 2853 500-6 2,257 8-315 26,120 
107 13890 9.520 59-51 334-2 1,496 5,492 17,220 
1-08 7-219 41-02 224-1 9941 3,637 11,380 
1-09 2087 5-622 2869 666-0 2-426 7,575 
1-10 «1279888 4488 2028 103-0 446-6 1,618 5,043 
4834 891727 66-08 231-4 708:3 
1-20 | 1259 | 1508 2-127 | 4-322 12-23 87-73 110-2 
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Carbon tetrachloride /Acetic acid azeotrope 


The effect of pressure on the composition of the azeotrope 


C. O. Havcurox 
Courtaulds Ltd, Research Laboratories, Coventry 


Summary A method is described for determining the composition of the carbon tetrachloride / 
acetic acid azeotrope at pressures within the range 10-1,000 mm Hg absolute. Reduction in 
pressure descreases the percentage of acetic acid in the azeotrope. The composition at 760 mm 
absolute is found to be 1-54% acetic acid by weight which is lower than previously reported 
values. 


Résumé—L auteur décrit une méthode pour déterminer la composition de l’azéotrope du 
mélange acide acétique — tétrachlorure de carbone par une échelle des pressions de 10 & 1000 mm 
: de mercure. L'abaissement de la pression diminue le pourcentage d’acide acétique dans 
55 lazéotrope. La composition 4 760 mm absolus est de 1,54°% en poids d'acide acétique, plus 
basse que les valeurs précédemment indiquées. 


INTRODUCTION 


The existence of an azeotrope of carbon 
tetrachloride and acetic acid boiling at 76-5°C 1300 
under 760 mm Hg absolute pressure and contain- 
ing 3.0%, acetic acid by weight has been reported 
by M. Lecar [1] and is published in the list 
of azeotropes complied by L. H. Horsiey [2]. 

The work described in this article was carried 
out to determine the effect of pressure on the 
composition of the azeotrope over the range 
10-1,000 mm Hg absolute. 


EXPERIMENTAL METHOD 


Six mixtures A-F were made up containing 
400 g of pure carbon tetrachloride and approxi- 
mately 0.5, 2, 4, 6, 8 and 10 grams respectively 
of glacial acetic acid (99-8% by freezing point). 
These mixtures were distilled in turn from a 
500 ml. flask through a Dufton spiral glass 
fractionating column 9-0 inches high x 10 
inches diameter, still head, and condenser, the 
distillates being collected in a receiver. All 
connections were made with ground glass joints. 
The receiver had a side arm connected to a 
source of vacuum and a mercury manometer. 

Each distillation was carried out at atmospheric yg, 1, Effect of pressure on the composition of carbon 
pressure and at three or four different reduced tetrachloride—acetic acid azeotrope. 
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Table I. Acetic acid contents of distillates and residues from miatures A-F. 


Pressure BP. Mis Weight %, acetic acid in 
(mm) distilled 
MIXTURE A 761 76-0 12 ‘B10 107 
364 19 173 084 
211 39-0 15 19 078 
137 20 22 “084 O71 
112 24 18 “O79 073 
MIXTURE B 753 | 75-6 15 “798 5385 
538 65-0 24 “687 
347 52:5 37 565 “557 
209 38-8 29 573 
MIXTURE C 764 76-0 27 1-20 0-98 
559 66-0 4 1.02 101 
421 57-7 16 0-93 1-02 
288 47-1 24 0-83 1-05 
144 _ 12 0-65 1-02 
MIXTURE D 753 75-6 17 1-54 1:54 
614 | 69-2 22 1-41 1:56 
397 | 56-4 2 1-28 1-60 
MIXTURE E 753 75:8 15 201 2-39 
860 80-8 207 2-37 
MIXTURE F 753 75-8 0 | 2-20 3-10 
648 | 20 2-08 
491 62-6 20 1-92 8-12 
231 41-5 30 171 -- 
| 
pressures. About 20 mil. of distillate were Table II. Composition and boiling point of 
collected at each pressure and the distillate azeotrope. 
and residue were both analysed by titration of 
5 or 10 ml. with 0-1 N sodium hydroxide. The — BP. Weight % 
compositions of distillates and residues were (mm) ec) acetic acid 
plotted against pressure and the point where 
the two curves intersected was taken as the 1400 — 30 
composition of the azeotrope. 
530 64-6 1-00 
musvate 340 515 0-58 
The results of the six series of distillations 90 18-7 0-72 
are given in Table I and the distillate and residue 


compositions are plotted against pressure in the 
figure. A curve drawn through the points where 
the distillate-residue compositions intersect gives 
the composition of the azeotrope. 


In experiments E and F intersections of the 
distillate and residue composition lines lay 


outside the operating pressure range of the 
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upparatus and had therefore to be found by 
extrapolation, with consequent reduction in the 
accuracy of the results. 

The boiling points of the six azeotropes thus 
established have been obtained from the tempera- 
ture-pressure data in Table I by plotting boiling 
point against pressure for each of the mixtures 
A-F and reading from the curves the boiling 
point of the azeotrope at its corresponding 
pressure. Table 2 gives the composition and 
boiling point of the azeotrope thus found. 


CONCLUSIONS 
1. The composition of the carbon tetra- 


chloride-acetic acid azeotrope at atmospheric 
pressure is 1-54% acetic acid by weight. This 
disagrees with the previously reported figure 
of 3% by weight published by Lecar [1]. 

2. Reduced pressure decreases the percentage 
acid in the carbon tetrachloride /acetic acid. 
From the slope of the curve in Fig. 1 it would 
appear that an azeotrope is no longer formed 
at pressures below about 50 mm Hg absolute. 
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Book reviews 


H. Mresswer and U. Gricuct (Editors): Fortschritte 
der Verfahrenstechnik 1952/53. (Progress in Process 
Technique in 1952-53), Verlag Chemie, GMBH, Weinhem- 
Bergstr. 1954. 496 pages. Price DM 39. 


Tuts is the first volume of a proposed series, with succeeding 
volumes to be published by the Farbenfabriken Bayer at 
intervals of one to two years, in which an attempt will be 
made to review and summarize the literature concerned 
with unit operations and process techniques. In general 
each volume will be limited to the literature which has 
appeared since publication of the preceding volume, al- 
though in the first volume references are made to works 
published before the period of 1952-53. The book is 
divided into sections covering the usual unit operations, 
such as fluid flow, heat transfer, gas absorption, and dis- 
tillation; the techniques involved in the use of high 
temperatures, low temperatures, high pressures, and 
vacua; and reactors, materials of construction, and 
corrosion. Each section was written by one or more 
specialists who are well known in their respective fields. 
Very extensive bibliographies are presented, in which, 
somewhat unfortunately, the references are listed in order 
of discussion rather than alphabetically. Author and 
subject indexes are also provided. 

The book is well written and the various fields appear 
to have been thoroughly covered. The reviews are critical 
and are not simply a collection of abstracts. Many illus- 
trations of equipment and graphs presenting technical 
data are included. The book should be particularly useful 
to one who wants to keep up with progress in the unit 
operations in a general way or to one who is making a 
literature survey in a particular field. 

Since reviews of unit operations, unit processes, and 
materials of construction appear regularly in publications 
written in English, such as “ Industrial and Engineering 
Chemistry," Fortschritte der Verfahrenstechnik will probably 
be of limited interest outside of Germany. 

Raywonn W. Sovraworta 


M.S. Perers: Elementary Chemical Engineering. 
pp. 322 + ix. McGraw-Hill Book Co. Inc. $6. 


Tue declared object of this book is to present a unified 
picture of chemical engineering for persons without 
experience or training in the subject. This is an ambitious 
undertaking, and the author has dealt with a wide range 
of subjects starting with elementary physical chemistry 
and stoichiometry, following with fluid flow, heat transfer 
and the unit operations, and concluding with chapters 
on chemical technology, economics and plant design. 


It is difficult to see how so many subjects could receive 
adequate discussion in a text of 300 pages. In the book 
under review the theoretical treatment of the unit 
operations has been condensed by omitting those aspects 
of the subject which require more than an elementary 
treatment. The basic equations of heat and mass transfer 
are presented with little or no discussion of the theoretical 
principles involved, though their practical application 
is considered more fully and illustrated with worked 
examples and problems. In some instances this process 
of condensation has been carried to a point of absurdity, 
for example solvent extraction is dealt with in two pages 
and chemical kinetics receives only half a page. A single 
chapter of twenty pages is devoted to the description 
of all types of chemical plant and equipment, ranging 
from pipe fittings and pumps to evaporators, towers. 
filters and grinding plant. 


The book could be recommended as a useful revision 
text-book for students receiving a course of tuition and 
lectures on chemical engineering. However, the uninitiated 
reader would inevitably gain the impression that chemical 
engineering comprises a smattering of chemistry, some 
elementary mathematical exercises and formulae relating 
to the unit operations, and very little of an engineering 
nature. 

W. S. Norman 


